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ABSTRACT 


The  tip  relief  effect  on  the  performance  of  helicopters  operating  at  high 
subsonic  tip  Mach  numbers  is  calculated.  The  theoretical  method  based  on 
the  complementary  ving  model  for  tip  relief  is  applied.  The  method  of 
application  in  complex  blade  element,  computer  programs  is  described  and 
results  are  presented.  Similarly,  the  technique  and  results  from  a  simpli¬ 
fied  application  for  use  in  energy  type  prediction  methods  are  given. 

Since  tip  relief  is  closely  connected  with  the  problem  oi  ^  ..dicliiig  the 
compressible  power  loss  of  rotors,  the  techniques  for  predicting  this  power 
loss  are  reviewed  and  discussed  in  detail. 
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SUBSCRIPTS 


INTRODUCTION 


Compressibility  effects  have  long  been  recognized  a3  having  an  important 
influence  on  the  performance  of  high-speed  helicopters.  During  the  1960's, 
it  gradually  became  recognized  that  one  could  "live  with"  compressibility 
effects  on  rotors.  Earlier  it  van  found,  first  on  propellers  and  then  on 
rotors,  that  adverse  effects  of  compressibility  occurred  at  higher  Mach 
numbers  thar  expected  from  fixed-wing  results.  The  next  obvious  Btep 
was  to  attem.pt  to  alleviete  compressibility  effects  through  such  techniques 
as  the  use  of  thinner  01  supercritical  airfoil3  near  the  tip  of  the  blade, 
snaping  of  the  tip,  etc-  Parallel  efforts  were  directed  at  developing 
techniques  for  predicting  performance  of  rotors  considering  compressibility 
and  correlating  flight  test  data.  Further  studies  and  experimental  ior. 
are  ..or.  tinning  at  a  fast  psce. 


This  report  presents  the  results  of  ar  mve  tigation  conducted  by 
Aerophysics  Company  to  obtain  a  better  theoretical  understanding  of  com¬ 
pressibility  effects  or,  helicopter  rotor  performance.  It  exploits  earlier 
findings  by  leNard*  based  or.  a  hypothesis  suggested  by  Trant  tnat  a  theo¬ 
retical  potential  flow  method  of  source  -  sink  cingularities  similar  to 
the  method  developed  by  Anderson  for  fixed  wings2  could  account  for  the 
effect  of  the  compressible  three-dimensional  relief  on  the  torque  required 
for  a  helicopter  rotor.  Results  of  Reference  1  were  systematically 
analyzed  in  the  study  reported  herein  for  application  with  existing  rotor 
performance  analyses.  Cince  the  main  point  of  the  study  i3  tip  relief, 
i.e.,  a  relief  effect  net  accounted  for  in  conventionally -used  compressi¬ 
bility  calculations,  it  was  nececsary  to  carefully  review  these  conven¬ 
tional  compressibility  analyses.  Those  analyses  are  summarized  in  the 
first  part  of  the  report.  Then,  a  brief  review  cf  the  tip  relief  analysis 

of  Reference  1  is  presented.  This  tip  relief  analysis  is  then  applied  and 
used  wizh  existing  blade  slement  computer  programs,  both  hover  and  forward 
flight,  and  in  a  very  simplified  biaie  element  analysis  whereby  the  torque 
decrease  due  to  tip  relief  is  calculated  as  a  separate  <  ^mponent  of  the 
total  torque.  Numerics i  examples  are  calculated  and  performance  charts 
of  the  tip  relief  effect,  usable  for  design  purposes,  are  drawn. 


As  some  of  today's  helicopter  reters  sne cur.t; 
transonic  similarity  rules  derived  for  fixed  wings  were  examined.  It  was 
hoped  that  these  similarity  lavs  could  be  extended  to  rotors  and  combined 
with  the  tip  relief  theory.  Although  this  subject  was  examined  and  is 
discussed  in  this  report,  it  has  not  yet  been  successfully  applied  to 
rotors . 


Exper imentai  verification  of  the  Anderson  -  LeNord  theory  for  tip  relief 
cf  rotors  ias  been  considered;  however,  no  simple  and  relatively  inexpen¬ 
sive  method  has  been  suggested  that  will  effectively  isolate  the  tip 
relief  effect  tc  the  precision  required  to  verify  the  theory. 

It  as  concluded  that  the  Anderson  -  I.eNard  tip  relief  method  has  been 
success!  .illy  reduced  to  practice.  Attention  should  now  be  concentrated 
on  the  other  aspects  of  helicopter  rotor  compressibility. 


COMPRESSIBILITY  EFFECTS  IN  HELICOPTER  AERODYNAMICS 


The  problems  of  high-speed  rotors  and  propellers  have  been  under  investi¬ 
gation  since  the  early  19^0's.  Propeller  compressibility  problems  became 
acute  on  high-speed  fighters  during  World  War  II.  A  number  of  approaches 
to  subsonic  and  supersonic  propeller  performance  evolved.  During  this 
period,  experiments  showed  a  discrepancy  between  the  test  results  and 
blade  element  theory,  in  that  the  tests  gave  an  apparently  lower  drag 
divergence  Mach  number.  Many  investigators  have  attempted  to  find  means 
of  predicting  the  compressible  flew  effects  as  are  discussed  here. 

At  high  forward  speeds  it  is  difficult  to  separate  compressibility  and 
stall  effects,  as  both  often  occur  at  the  same  time  on  the  rotor.  At 
the  advancing  tip,  the  forward  speed  and  the  rotational  speeds  add, 
resulting  in  speeds  approaching  Mach  one.  But  on  the  retreating  side  the 
two  speeds  subtract,  resulting  in  a  region  of  reverse  flow  and  areas 
where  the  flapping  blades  are  at  or  above  the  stall  angle  of  attack. 

The  dividing  line  between  stall  losses  and  compressibility  losses  is 
seldom  deducible  from  test  data,  and  both  areas  must  be  taken  into  account 
in  rotor  analysis. 

COMPRESSIBLE  FLUID  DYNAMICS  OF  ROTARY  WINGS 


The  motion  of  rotating  wings  through  a  fluid  should  be  describable  by 
Cue  Ltttoic  differential  equations  ol'  fluid  metier..  But  the  complexity  of 
the  rotary-wing  problem  is  such  that  only  recently  has  any  real  progress 
been  made  in  the  basic  theory.  The  fluid  dynamic  equations  governing  the 
motion  of  the  fluid  for  a  rotating  wing  can  be  derived  from  the  equations 
of  continuity,  momentum  and  energy.  This  derivation  is  presented  in  a 
number  of  places . 1 » 3 >4 *  5 *6 


Sears7  made  or.e  of  the  first  attempts  at  beginning  rotary -wing  theory 
with  the  basic  governing  equations  of  fluid  motion  in  1950.  He  investi¬ 
gated  the  flow  potential  function  of  an  infinitely  long  cylindrical  blade, 
rotating  in  incompressible  flow.  The  result  indicated  that  the  potential  per 
unit  velocity  in  the  plane  transverse  to  r.bp  rotating  cylinder  is  the 
same  as  for  a  similar  cylinder  in  steady  plane  flow.  Further,  immediately 
adjacent  to  the  cylinder  (at  least  for  the  circular  cylinder  cane  given 
in  Sears'  paper)  the  flow  is  only  in  this  transverse  plane.  The  con¬ 
clusion  which  may  be  drawn  is  that  the  pressure  coefficient  and  therefore 
the  force  coefficients  from  potential,  incompressible  flow  are  the  same 
aa  that  from  blade  element  theory.  Thus  the  validity  of  t be  blade  element 
approach  was  indicated,  and  it  has  been  used  for  performan  calculations 
and  theoretical  work  including  turbulent,  compressible  rotor  boundary 
layers  (for  example,  see  hickr  end  Nash,  or  Clark  and  ArnoldJ,9).  Only 
in  December  1970  was  it  shoVn  by  Goorjian  and  McCroskey10  that  tne  Sears 
approach  is  not  necessarily  valid  when  the  rotating  blade  is  of  finite 
span.  The  major  theoretical  effort  over  the  past  15  years  has  been 
toward  determining  the  effect  of  the  trailing  vortices.  The  trailing 
vertices  give  the  lue&I  Inflow  velocity,  wmch  Is  a  major  quantity  in 
determining  the  local  angle  of  attack,  particularly  near  the  tip  and  in 
hover.  A  great  amount  of  literature  is  available  cn  this  approach* 1 >12 


2 


and  thus  will  not  be  discussed  in  detail  here,  but  it  generally  considers 
incompressible  flows. 

In  the  case  of  uniform  flows,  the  well-known  Prandtl-Glauert  transformation 
may  be  used  to  change  the  equation  for  small  pertubation  velocities  in  com¬ 
pressible  flows  to  the  equation  for  ar.  equivalent  incompressible  flow.  It 
was  shown  in  Reference  5  that  small  perturbation  equations  for  a  rotating 
blade  do  not  reduce  to  a  simple  form,  amenable  to  the  Prandtl-Glauert 
transformation.  Further  assumptions  are  needed  so  that  a  Prandtl-Glauert 
type  of  transformation  may  be  used,  as  shown  in  Reference  1.  The  differ¬ 
ence  between  the  two  cases  is  that  for  uniform  flow  the  compressible, 
small-pertubation  equation  has  constant  coefficients,  while  in  rotating 
flow  it  is  again  a  linear,  second-order  partial  differential  equation, 
but  the  coefficients,  are  now  functions  of  the  space  variabl.es. 

There  have  been  two  recent  papers  concerning  the  compressible  potential 
flow  over  nonlifti ng ,  hovering  rotor  blades.  The  earlier  one  by  Sopher5 
derives  the  equation  for  a  source  distribution  representing  a  thin  airfoil, 
rotating  blade  of  finite  span  in  the  subsonic,  compressible  regime.  The 
ir.visoid  velocity  distribution  was  then  calculated  at  several  spanwise 
stations.  Results  from  numerical  solutions  to  the  equations  for  transunic 
flow  over  a  rotating  blade  were  presented  by  Caradonna  and  Isom.6  This 
paper  indicated  that  as  the  Mach  number  increases,  "the  nonlinear  behavior 
of  the  wing  near  the  tip  becomes  more  severe  than  that  of  the  rotor". 

Thus  linearized  theory  (subsonic  flow)  has  validity  t.o  higher  subsonic 
Mach  cumbers  for  a  rotor  than  for  a  wing. 

The  asymptotic  expansion  to  the  linearized  subsonic  potential  equation 
in  terms  of  the  distance  from  the  tip  was  derived  by  Caradonna  and  Isom.6 
To  a  first  approximation  near  the  tip,  the  flow  was  found  to  be  the  same 
as  the  three-dimensional  flow  in  the  tip  region  of  a  semi-infinite  wing 
moving  at  a  uniform  speed  equal  to  the  blade  tip  Mach  number.  This  is 
valid  for  high-aspect-ratio  rotors  where  the  Mach  number  gradient  is 
large.  But  in  the  high-subsonic  compressible  range,  small  enunges  in  the 
Mach  number  car,  result  in  large  changes  in  the  drag.  Thus  this  approxi¬ 
mation  must  be  used  with  care. 

EARLY  F.ypERIMB'MTAL  WORK 

Although  seme  experimental  work  was  done  during  the  19^0's,  by  Gutsche13 
and  Gustafson14  among  others,  the  major  work  in  this  area  is  the  result 
of  a  series  of  tests  conducted  on  the  NACA  Langley  test  tower  under  the 
direction  of  Paul  Carpenter. 

The  Carpenter  tests  were  performed  on  a  number  of  different  rotors, 
varying  blade  airfoil  section.  Reports  of  these  tests  are  given  in 
references  15  through  25.  Comparison  cf  the  test  results1 6 t20  with  strip 
analysis  of  Reference  ?6  using  two-dimensional  airfoil  data  indicated 
that  the  theory  overestimated  the  compressibility  drag  rise.  Carpenter21 
gave  as  reasons  for  these  discrepancies  the  effect  of  centrifugal 
pumping  and  the  drag-alleviating  effects  of  three-dimensional  flow  at 
the  blade  tin  ting  from  c.  reduction  uf  local  Mach  number.  This  led 

him  to  formulate  the  concept  of  synthesized  blade  section  data. 
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Synthesized  data  is  obtained  by  adjusting  the  airfoil  Bection  properties 
until  the  strip  analysis  performance  prediction  agrees  with  the  measured 
performance.  This  is  a  trial  and  error  approach,  but  certain  guidelines 
from  theory  and  test  axe  used  to  obtain  reasonable  values  of  the  airfoil 
section  properties.21 

The  synthesized  data  is  obtained  from  hover  teats  on  a  test  tower. 
Reasonably  large  blades  (26.8  feet  radius  and  16.^  inches  chord  for  the 
NACA  0012  airfoil,  for  example)  were  used,  to  obtain  high  enough  Reynolds 
numbers,  but  surface  condition  differences  may  be  important.  Concerning 
application  to  forward  flight,  Carpenter  Btates  that  the  synthesized 
data  should  be  valid  until  proved  otherwise.21  Tc  date  this  hypothesis 
has  been  neither  proved  nor  disproved. 

One  of  the  possible  applications  envisioned  for  the  tip  relief  theory 
is  to  provide  a  link  between  synthesized  airfoil  data  and  two-dimensional 
airfoil  data.  The  synthesized  airfoil  data  is,  in  a  sense,  an  average 
of  the  two-dimensional  data  with  tip  relief  correction. 

BLADE  ELEMENT  ANALYSIS 


Blade  element  or  strip  analysis  forms  the  basis  of  almost  all  performance 
and  rotor  dynamics  calculations .  The  so-called  ’’energy"  methods  are  in 
most  instances  the  result  of  simplifying  the  strip  analysis  so  that  the 
integration  may  he  carried  out  analytically  (i.e.,  resulting  in  an 
algebraic  equation),  rather  than  requiring  numerical  computation. 

The  basis  of  blade  element  theory  is  the  assumption  that  each  element 
can  be  considered  aerodynamic ally  as  an  independent  two-dimensional 
airfoil  segment.  The  local  forces  are  calculated  from  the  resultant 
velocities  at  the  element,  and  the  total  for"®*  are  found  by  integrating 
radially  and  azimuthally.  From  the  forces,  the  motion  of  the  element  is 
determined.  But  the  motion  (especially  flapping)  has  an  effect  on  the 
local  velocities  and  therefore  the  forces,  requiring  iterations. 

The  early  performance  methods  made  a  number  of  assumptions  to  avoid  the 
iterative  calculations.  These  concerned  the  induced  velocities,  the 
lift  and  drag  coefficients,  the  flapping  motion,  etc.  An  example  may  he 
found  in  Chapter  8  of  the  hook  by  Gessow  and  Myers,27  although  other 
workers  have  used  different  sets  of  assumptions. 

With  the  development  of  computers,’  an  iterative  method  is  used,  which 
in  essence  was  first  outlined  by  Gessow.26  There  have  been  a  number  of 
improvements  made  in  various  areas  since  that  time.  The  most  widely 
used  results  of  this  technique  are  the  tables  and  charts  prepared  by 
Tanner.2 “»29 

There  is  much  work  still  being  done  and  controversy  over  the  airfoil  data 
to  be  used  for  these  computerized  strip  analyses.  The  latest  suggestions 
concern  the  problem  of  unsteady  stall  hysteresis.30  The  Tanner  charts26 
use  airfoil  data  from  "two-dimensional,  wind  tunnel  tests  of  a  production 
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blade  specimen  using  the  RACA  0012  airfoil  section".  A  similar  set  of 
charts  was  prepared  using  the  Carpenter21  synthesized  data  by  Kisielovski 


et  al 
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To  our  knowledge  there  has  been  no  meaningful  comparison  of 


these  two  prediction  methods,  although  both  appear  to  be  widely  used. 


The  method  proposed  to  account  for  tip  relief  is  to  use  the  two-dimensional 
airfoil  data,  but  to  make  the  corrections  as  indicated  by  the  Anderson  - 
LeHard  theory.  Detailed  derivation  of  the  tip  relief  theory  is  given  in 
Reference  1,  while  the  technique  in  general  and  the  results  of  use  in 
strip  analysis  are  described  in  the  following  chapter. 

To  obtain  an  indication  of  the  compressible  power  loss,  Gessov  xnd  Crim32, 
using  compressible  airfoil  data  in  the  method  described  in  Peference  26, 
calculated  the  power  required  for  several  rotor  configurations .  The 
compressible  profile  power  loss  of  the  various  configurations  was  com¬ 
pared  ct  equul  values  of  the  increment  in  tip  Mach  number  above  the  drag 
divergence  Mach  number.  The  results  bhoved  that  the  different  config¬ 
urations  had  approximately  the  same  loss  in  profile  power  coefficient  to 
solidity  ratio  when  compared  in  this  manner.  This  suggested  that  the 
increment  in  power  coefficient  to  solidity  ratio  due  to  compressibility 
depends  only  on  the  difference  betveen  the  drag  divergence  and  the 
advancing  tip  Mach  numbers.  A  more  detailed  discussion  is  given  later 
in  this  report  in  reference  to  Figure  9* 

jju-iHUJAl  METHODS 


To  reduce  the  data  and  decrease  the  amount  of  flight  test  required,  per¬ 
sonnel  at  the  Edvards  AFB  Flight  Test  Center  have  been  developing  an 
approach,  mainly  empirical,  to  account  lor  compressibility  effectc. 

Drawing  on  fixed-wing  experience,  the  power  required  is  divided  into  two 
portions:  (l)  the  incompressible  power  and  (2)  an  increment  in  power 

due  to  compreBsible  flow,  which  depends  on  the  difference  between  the 
critical  Mach  number  of  the  blade  section  and  the  Mach  number  of  the 
element.  The  correlation  is  established  empirically  using  UH-1F,  CH-UyA 
and  CH-3E  flight  test  data.-3*3**  References  35,  36,  and  37  give  detailed 
flight  test  data  for  these  helicopters  flying  in  the  compressible  riow 
regime.  later  reports38*39  include  the  effects  of  power  losses  due  to 
retreating  blade  stall. 

The  early  empirical  work33*39  attempted  to  relate  the  pover  increase  due 
to  compressibility  to  the  value  of  the  advancing  tip  Mach  number  above 
the  critical  Mach  number.  These  works  differ  from  Reference  32,  where 
the  drag  divergence  Mach  number  is  UBed  instead,  Good  correlation  using 
test  data  from  two  different  helicopters  has  been  reported.31*  The 
resulting  empirical  relation  for  the  compressible  pover  rise  has  been 
used  with  varying  success  at  Edvards  and  elsevhere. 

More  recently,  a  nuaber  of  analytical/semi empirical  techniques  were 
developed  and  evaluated.30*-9  These  techniques  used  various  approximate innn 
uu  the  drag  in  compressible  flow  and  the  lift  in  the  stall  region.  Thus 
both  the  compressibility  and  the  stall  regions  are  considered.  The  basic 
theory  was  modified  empirically  to  obtain  better  correlation  with  data. 
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The  techniques  consider  a  numbei  of  variables,  including  lift  coefficient 
at  the  advancing  and  retreating  tip,  lift  curve  slope  as  a  functioi  of 
Mach  number,  an  empirical  critical  Mach  number,  maximum  lift  coefficient, 
Mach  number  distribution  over  the  areas  affected  by  compressibility  and 
stall,  and  the  calculation  of  these  areas.  Results  are  compared  to  the 
CH-3C  flight  test  data  in  Reference  39.  Since  this  report  is  i.'idely 
available,  and  a  number  of  techniques  are  presented,  we  will  not  discuss 
it  in  detail.  Somsel3y  obtains  correlation  within  tlx  10“  4  in  power 
coeffic  ient  over  of  the  flight  e»  elope.  But  the  remaining  20$ 
occurs  at  the  higher  advance  ratios,  thrusts  and  Mach  numbers. 

ANALYTICAL  CALCULATION  OF  THE  COMPRESSIBILITY  EFFECTS  ON  ROTOR  PERFORMANCE 


B>  analytical  calculation  we  mean  those  approaches  which  attempt  to  yield 
an  algebraic  expression  for  the  forces  on  a  rotor  due  to  compressibility 
and  thus  separate  the  forces  due  to  compressibility  from  the  other  com¬ 
ponents  of  the  total  force.  This  falls  into  the  realm  of  "energy"  type 
approaches.  These  analytical  calculations  use  one  of  two  methods: 
actuator  disc  theory  or  blade  element  analysis. 

In  the  late  19^0's  and  early  ISiO's,  the  problem  of  high-speed  propellers 
vaa  investigated  using  actuator  disc  theory  in  compressible  flow . 4^ »4 1 »42 
In  these  investigations,  the  compressible  flow  through  a  contracting 
stream  tube  with  energy  added  at  the  actuator  disc  is  considered.  This 
flow  model  represents  axial  flight  and  is  more  applicable  to  axial  turbines 
and  "o  propellers  at  high  forward  speeds.  The  problem  of  compressible 
tii  speed  doe6  not  appear  in  the  results  due  to  the  infinite  number  of 
blaae  assumptions  inherent  in  the  actuator  disc  approach.  This  actuator 
disc  approach  was  combined  with  blade  element  theory  by  Laitone4j  for 
application  to  helicopters. 


An  approach  which  appears  to  be  more  applicable  and  useful  was  developed 
by  Head,44  by  the  use  of  theoretical  and  empirical  rules  for  the  variation 
of  the  drag  and  lift  coefficients  with  Mach  number.  These  variations  are 
used  in  the  blade  element  equations,  which  are  then  integrated  to  find 
the  rotor  thrust  and  torque  coefficient  ratios  (i.e.,  the  ratio  of  com¬ 
pressible  to  incompressible  coefficients).  Further,  Head  considers  the 
Increase  due  to  compressibility  in  induced  power,  which  is  not  done  elsewhere. 


Head's44  paper  may  be  outlined  as  follows:  First,  Head  hypothesises  a 
variation  of  lift  coefficient  with  Mach  number.  This  variation  was  then 
used  in  the  equations  for  torque  and  thrust  for  the  hover  case  ,  which 
vaa  then  solved.  He  then  examined  the  forward  flight  case.  The  com¬ 
pressible  profile  drag  power,  for  both  hover  and  forward  flight,  was 
then  considered,  hypothesizing  a  variation  of  drag  coefficient  with  Mach 
number.  This  outline  is  followed  In  the  discussion  below. 


The  lift  coefficient  is  assumed  by  Head  to  vary  in  accordance  with  the 
Prandtl-Glauert  law  for  Mach  numbers  up  to  lift  divergence.  Above  lift 
divergence,  a  linear  decrease  with  Mach  nunber  is  assumed. 
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Making  the  usual  assumptions ,  such  as  small  inflow  angles,  neglecting  the 
effects  of  drag  on  the  thrust  and  induced  torque,  etc.,  the  compressible 
to  incompressible  ratio  of  thrust  and  induced  torque  eocfficic  nts  sure 
calculated  for  the  hover  cose.  These  ratios  are  functions  of  e  tip 
angle  of  attack,  lift  divergence  and  rotational  Mach  numbers,  and  the 
variation  of  lift  coefficient  with  angle  of  attack  and  Mach  number. 

For  the  hover  case,  it  is  found  that  not  only  the  torque  but  the  thrust 
increases  with  Mach  number.  Thus  most  of  the  increased  induced  torque 
due  to  compressibility  is  due  to  the  increase  in  thrust  due  to  com¬ 
pressibility.  Thus,  given  a  required  thiust,  one  would  reduce  the  angle 
of  attack  (actually  collective  pitch)  to  make  Cpc./CTinc  equal  to  unity. 
With  this  condition,  the  increase  in  induced  torque  due  to  compressibility 
does  not  exceed  3?  at  and  decreases  above  the  lift  divergence  Mach 
number,  since  there  is  a  relieving  effect  <-  \  lift  that  results  in  the 
inboard  movement  of  the  centroid  of  the  lewd.  These  results  from 
Reference  1*1  are  shown  in  Figure  1.  Part3  (a)  and  (b)  of  Figure  1  slow 
the  effect  of  compressibility  on  thrust  and  argle  of  attack.  In  part  (c), 
it  can  be  seen  that  the  induced  torque  due  to  compressibility  at  constant 
thrust  is  very  snail  and  nay  be  neglected,  while  at  a  constant  angle  oi 
attack  this  is  not  the  case. 


These  results  indicate  that  in  hover,  calculating  only  the  incompressible 
values  of  the  thrust  and  induced  torque  is  a  reasonable  assumption  in 
the  usual  energy  calculations.  The  major  effect  of  compressible  flow  on 
rotor  performance  is  upon  the  profile  power  a  id  upon  the  airfoil  angle  of 
attack.  Much  of  the  other  work,  as  described  herein,  has  been  concerned 
with  the  profile  power  due  to  compressibility.  On  the  other  hand,  the 
effect  of  compressibility  on  the  angle  of  attack  has  been  neglected. 

For  forward  flight..  Head  states  that  the  effect  cf  compressibility  on  the 
induced  power  nay  be  neglected  and  that  the  equation  for  thrust  cannot 
be  analytically  integrated  but  must  be  numerically  integrated,  which  he 
iues  not  do. 


The  calculation  of  compressible  profile  drag  power  is  based  on  the 
assumption  that  the  section  drag  coefficient  variation  vith  Mach  number 
may  be  represented  by  a  parabola.  Head  writes 

AC  ■  /  /  (w  Cob  +  r)3  AC^  dr  (3) 


where  the  limitr  of  integration  on  and  r  are  chosen  to  represent  the 
"portion  of  the  rotor  disc  operating  above  the  drag  divergence  Mach 
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Figure  1.  Effects  of  Compressibility  on  Rotors  in  Hover. 


number"  The  incremental  drag,  ACn,  is  the  increase  in  profile  drag  due 
to  compressibility.  Head  assumestiasically  a  parabolic  variation  of  ACp 
with  Mach  number,  but  makeo  further  correction,  auch  a a  an  increase 
of  0.06  for  the  rotor  drag  divergence  Mach  number  above  the  two-dime ns ional 
value,  quoting  Reference  l1*.  He  carries  out  the  integration  indicated 
in  equation  (3)  for  both  hover  and  forward  flight.  The  results  are  pre¬ 
sented  in  chart  form,  and  ACn  ,  may  be  calculated  given  Mp  jq  p,  M^d  and 
the  parameters  for  the  parabolic  variation  of  AC^  vs.  M. 

Results  for  hover  appear  very  good  vhen  the  Carpenter15  data  is  compared  to 
tiie  theory.**1*  In  forward  flight, ,  some  limited  comparison  made  to  the  data 
of  Reference  36  appears  to  give  low  values  of  the  AC.^,  but  this  result 
should  not  be  taken  as  conclusive. 

Further  work  with  this  theory  is  necessary  to  determine  the  importance  of 
tiie  change  in  angle  of  attack  due  to  compressibility ,  the  change  in  induced 
torque  in  forward  flight  due  to  compressibility ,  and  comparisons  to  test 
data. 
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Til'  RELIEF  IN  ITT-  RKLATluRElilE  TO  COMER EEETJ31L1 TV 


The  existence  of  a  tip  relief  effect  has  beer,  suggested  by  various 
sources,  as  discussed  in  the  previous  chapter.  The  tip  relief  effect  ia 
o  difference  between  the  real  tluee-dimcr.sional  flow  over  a  finite  spun 
wing  or  rotor  blade  and  the  calculation  for  this  wing  or  rotor  blade  nude 
using  two-dimensional  airfoil  data.  In  this  chapter  the  theory  arid 
lesults  fi  cm  one  approach  to  calculating  the  magnitude  of  the  tip  relief 
are  given.  The  equations  of  the  theory  were  derived  in  detail  in 
Reference  1.  Here  the  basic  elements  of  the  theory  for  its  use  are 
given  and  the  results  from  its  application  in  several  performance  cal¬ 
culations  . 

DESCRIPTION  OK  THEORY 


For  fixed-wing  aerodynamics,  the  major  effect  of  aspect  ratio  occurs  in 
the  spunwiee  variation  of  induced  velocity  and  the  resulting  induced  drag 
and  angle  of  attack.  This  is  all  due  to  the  production  of  lift  by  the 
finite  fixed  wing.  At  lower  Mach  numbers  this  approach  is  satisfactory, 
but  at  higher  subsonic  Much  numbers  a  further  correction  is  needed,  as 
was  described  by  Anderson." 

VI. on  u  tnree-dimenaionul  object,  such  as  a  wing,  is  pluced  in  a  uniform 
flow,  the  stream  lines  around  the  body  may  deviate  from  the  free  stream 
in  both  directions  perpendicular  to  the  wing  chord  line.  For  a  iin.il ar 
two — dimensional  wing,  o*.Tj  ucviuiiun  in  the  uiieeciun  perpenu lculur  lo 
the  plane  Including  both  the  chord  and  the  span  is  possible.  Thus  one 
can  qualitatively  conclude  that  strcamwisc  disturbances  produced  by 
tlu'ce-dimensional  bodies  are  lever  in  magnitude*  than  those  produced  by 
two-dimensional  wings.  In  the  case  of  a  r«  ct-angulur  pi  an  form  wing  of 
hish  uni ect  ratio,  at  the  midspan  the  disturbances  are  nearly  the  some  as 
for  two-dimensional  flow  and  diminish  us  the  tip  is  approached,  hence  the 
name  tip  relief. 

This  qualitative*  concept  was  quantified  for  fixed  wings  by  Anderson" 
through  the  use  of  linearized  wing  theory  and  the  complementary  wing 
concept.  lit  showed  that  the  difference  between  the  two-  and  three- 
dimensional  characteristic  of  a  wing  can,  by  proper  averaging,  be  repre¬ 
sented  by  a  change  from  the  actual  free-streum  velocity  to  an  effective 
free-stream  velocity.  This  will  result  in  a  change  (decrease)  of  the 
effective  Mach  number  and  the  dynamic  pressure,  hence  in  the  drag 
coefficient. 

I.i near i zed  wing  theory  represents  the  thickness  by  a  source-sink 
distribution  and  the  camber  and  angle  of  attack  by  a  vortex  distribution. 
The  former  gives  the  so-called  r.onlifting  case,  while  the  vortex  distri¬ 
bution  gives  the  lift.  In  a  qualitative  sense,  we  have  shown  that  tip- 
relief  Is  due  to  displacement i  hence  we  need  eonsidei  only  the  thickness 
as  given  by  the  source-sink  distribution. 
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The  velocity  induced  due  to  thickness  of  a  finite  wing  in  the  integral, 
over  the  pi an form,  of  a  cource-Glnk  dt strilmtion  determined  by  the  local 
chordviae  airfoil  slope.  The  potential  function  of  this  wing  is 
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This  same-  potential  function 
vill  result  if  this  finite 
span  wing  is  represented  by  a 
wing  of  infinite  span  leas 
"complementary  wings". 

The  complementary  wing  is  u 
ciiiilar  wing  extending  from 
the  tips  of  the  finite  wing 
to  infinity  on  both  sides. 

In  the  equation  above  the 
limits  of  the  integral, 
using  the  coordinate  system 
of  Figur  e  2,  is  +e  and  -c 
on  and  -2A  and  0  on  n. 

In  r tyrn i >  1  t »r» »■> 4  »»_»« y  vili£ 

re*preaentation  the  potential 
is 


Figure  2.  Complementary  Wing  and  Coordinate 
System  for  a  Finite  King. 
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The  first  integral  is  identically  the  potential  oj'  an  infinite  or  two- 
dimensional  airfoil,  and  car.  be  represented  by  the  symbol  $x  s<*>  ■  The 
last  two  integrals  give  values  which  vary  with  position  tluoughout  the 
flew  field.  Let  -  G  U  ,x ,y  ,  l  }  /li.  represent  the  third  and  fourth  terms 
of  equation  (6).  Thus 
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Thee 


i  -  G(X,x,y,x) 


or  in  terms  of  the  local  velocities 


u  * 


2l 

3x 


U  ♦ 


<R 


u,  -  U  G  (X,x,y,z) 
X-CC  “X 


(8) 


■  Ujl-CMA.x.y,*})  +  ^  x»o> 


If  an  uverage  v<ilue  ol  can  be  found,  then 


u  »  U  (i-G  )  + 

®  x  X*«® 


(9) 


(10) 


and  the  velocities  (and  therefore  the  forces)  on  a  finite  wins:  can  be 
compared  to  the  velocities  a.id  forces  on  an  infinite  wing.  The  forces 
will  correspond  if 

U  (l-G  )  L,  ..  ■  U  I.  ..  ill) 

“  x  |  finite  “■  |  infinite 

wing  wing 

Thus  the  quantity  U®  Gx  represents  a  difference  or  change  in  tree-stream 
velocity  between  the  finite  and  the  infinite  wing  free-sti earn  velocities. 
With  pig'is  choser.  as  shewn. 


hU«. 

uT~ 


U 


CO 


| finite 
wing 


-  U 

oo 


infinite 

wing 


/  U 


(12) 


Anderson  in  Reference  2  did  carry  out  t)ie  averaging  process  over  the 
finite  wing  planform  to  obtain  GX(A)  for  use  in  fixed-ving  calculations , 
Because  of  the  spanvise  variation  of  flow  over  rotor  blades,  an  average 
which  is  aisc  a  function  of  spanwise  position  (or  distance  from  the  tip) 
would  be  more  useful  for  rotary  viiigs.  Therefore,  the  averaging  will 
be  carried  out  only  over  the  chord.  Since  thin  airfoil  assumptions  have 
been  made,  the  evaluation  may  be  done  in  the  ’*0  plane. 

_  +c 

G  (A,y)  =  (l/2c)  /  G^iX.x.y  ,zao)dx  (13) 

Then  the  two-dimensional  airfoil  data  now  used  in  blade  element,  analysis 
can  be  corrected  at  each  element  to  give  the  equivalent  tliree-dimensional 
data. 

Derivation  of  the  equations  and  the  mathematical,  manipulations  indicated 
are  shown  in  detail  in  Reference  1.  The  indiceted  integration  cannot  be 
done  anaiy t in  ally ,  and  a  Taylor  series  expansion  of  three  terms  is  used. 
For  helicopter  rotor  blades  of  the  usual  high  aspect  ratio,  orJ.y  the 
complementary  ving  near  the  tip  need  be  considered.  The  result  is 
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G  ^  iX*  =  r-  [  h  { c/y J  I  /c  +  hAc/y]  I  /c?  +  h  {c/y}  I  /c5  ]  (lU) 

X  U  oo  vTT  C  X  X  &  £•  j  J 

where 


h  {c/Y)  -  2  ’''(cATTU  — 1 

1  /[TTyT7-.! 


(15) 


h.lc/y)  -  i7  [U2  SkAhli?] 

-  l (c/y )2  +  l]5'//: 

h^tc/y}  = 


(16) 

(17) 


1_  [1,5.6  jjfil (°/y )S  +  (c/y)6]  1  63 (c/y) 14  +,  16< c/y)2  <  _8 
51  [(c/y)2  +  I]9/2 


T.ie  values  cf  In/t  do  not  v.j-y  appreciably  between  differing  airfoils,  ao 
they  ere  geometric  relatione  between  the  airfoil  profile  and  the 
rectangle  enclosing  the  profile.  Ine3e  values  calculated  for  a  parabolic 
arc  airfoiJ  arc 


I./t  -  -l.?33 
IJr  =  -C.600 
I3/T  =  -0.571 


This  change  in  velocaty  has  beer,  derived  for  incompressible  flows.  The 
Pranit.l-Giauert  transformation  is  applied  (in  the  manner  described  in 
pages  51  and  52  of  Reference  l),  result! iig  in  the  form 


AU  /U 
c 


00 


hn(6~,y} 


Vc 


2n-l 


(19) 


The  calculation  of  the  value  of  AUc/Uo,  has  been  programmed  into  a  simple 
subroutine  CORF  AC  described  in  the  appendix.  The  number  of  terms  of  the 
Taylor  series  used  car  be  varied.  There  is  only  a  slight  difference 
between  using  one  or  three  terms,  but  since  the  computer  time  is  only 
elightJy  increased,  it  is  suggested  that  three  terms  be  used,  by  adding 
the  result  Tumid  iron  replacing  the  variable  c/y  by  c/(2Ac-y) ,  since  2Ac 
defines  the  span  in  feet,  the  effect  of  the  opposite  complementary  wing 
is  found,  and  a  finite  rather  than  a  semi-infinite  wing  is  represented. 
The  use  of  a  finite  wing  does  not  appear  to  have  a  significant  effect.; 
thus  the  use  of  the  semi-infinite  representation  is  suggested. 

The  velocity  difference  can  be  expressed  in  terms  of  a  change  ir.  Kach 
number  and  in  -.rag  coefficient,  by  considering  the  compressible  quasi- 
one-dlmensional  flow  through  a  stream  tube. 
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(20) 


AM/M  =  (1  +  M2)  AU/U  =  r  AU  /U 

qd  >  oo  c  00 


acd/cd  =  - 


2-M 


1  +  M2 


AM/M 


(21) 


Thus,  given  the  two-dimensional  drag  coefficient  vs.  Mach  number  plot, 
the  three-dimensional  characteristic  can  be  found  by  using  the  above 
two  equations  as  sketched  in  Figure  3  (a).  We  should  note  that  the 
change  in  Mach  number  giver,  by  equation  (20)  causes  a  much  larger  change 
in  drag  coefficient  in  the  subsonic  drag  rise  range  than  the  AC^ 
correction  given  in  equation  (21). 

The  technique  as  derived  here  gives  the  difference  in  flow  character¬ 
istics  between  two-  and  three-dimensional,  inviscid  and  norlifting 
wing.-.  Since  it  is  the  change  that  is  calculated,  the  equations  can  be 
app-5 ed  to  viscoui  and  lifting  wings.  The  tip  relief  effect  is  pri¬ 
marily  a  di spire. o-tcnt  effect  which  is  represented  by  the  nonlifting  wing. 

V? !K  USE  07  THE  Til  RELIEF  CORRECTION  II.'  STRIP  ANALYSIS 


Strip  analysis  calculations  assume 
aerodynamically  as  an  independent 
velcci'tisc  Mslc^  ,  etc*  , 


that  each  element  can  be 
two-dimensional  airfoil  s 

4-  a  4*  4-  V  e.  T  r.  f\  4  a  *• 


considered 
egment  at 

^  ^  4-  U  ~ 

V/*  wuu 


element  on  the  rotor  disc.  The  tip  relief  correction  changes  the  Mach 
number  as  per  equation  (20),  and  the  drag  coefficient  is  then  found 
at  this  new  effective  Mach  number.  The  drag  coefficient  is  further 
corrected  by  the  ACjj  given  by  equation  (21).  No  change  is  made  for  the 
lift  or  moment  coefficients  since  the  major  difference  between  calcu¬ 
lations  and  test  is  in  terms  of  drag,  although  the  technique  can  be 
applied  to  these  other  coefficients. 


The  application  for  the  tip  relief  correction  can  best  be  described  by 
the  step-by-step  outline  given  below.  Instead  of  directly  determining 
the  drag  coefficient  as  is  now  done  in  strip  analysis,  the  following 
steps  should  be  carried  out: 


1.  Calculate  th  correction  factor  using  equation  (19)  for  the 
particular  b_ade  element. 

2.  Determine  M  =  M*,  +  AM  using  equation  (20).  (22) 


Determine  the  drag  coefficient  Cp.  but  at  M  as  is  row 

Deff,  M  eff 

done,  all  other  parameters  (a,  Ke,  etc.)  remaining  the  same. 

This  is  the  drag  coefficient  with  Mach  correction  only,  Meff 


1*.  Using  equation  (21),  calculate  the*  effective  drag  coefficient 


-•  c. 

‘el’.,  Meff 


( i  +  i,r  /r.  } 
h ■  v  ' 


(23) 


lb 


(a)  Determining  finite  wing  drag  variation  with  Mach 
number  from  infinite  wing  data.  Note  that  the 
process  may  be  reversed. 


Finite  Wing 

(b)  Determining  effective  of  finite  wing  (or  its 
element)  given  infinite  wing  data  and  Mach  number 
of  finite  wing  (or  its  element). 


Figure  3.  Application  of  the  Tip  Relief  Correction. 


The  Bigne  used  in  equations  (22)  and  (23)  need  sane  further  explanation, 
a s  they  nx\y  appear  inconsistent.  The  reason  io  that  the  derivation  of 
equations  (15)  through  (21)  is  made  to  go  from  two-dimensional  to  • 
three-dimensional  characteristics ,  or  vice  versa.  But  the  application 
described  here  is  that  given  the  three-dimensional  flight  parameters 
(i.e.,  Mach  number)  and  the  two-dimensional  aerodynamic  characteristic 
(namely,  drag  coefficient),  determine  the  three-dimensional  drag 
coefficient.  This  is  shown  in  Figure  3  (b). 

These  steps  have  been  incorporated  intfc  the  strip  analysis  computer 
programs  for  rotor  performance  calculation  that  are  used  by  the  Euatis 
Directorate,  USAAMRDL.  These  programs,  one  for  hover  and  another  for 
forward  flight,  use  straightforward  blade  element  strip  analysis 
techniques  (which  will  not  be  described  here).  These  programs  determine 
the  airfoil  characteristics  by  calling  an  external  subroutine  which  stores 
the  airfoil  data.  A  simple  option  was  added  to  the  main  programs  either 
to  call  the  airfoil  subroutine  as  vas  done  previously  or  to  calculate 
the  performance  using  tip  relief  by  calling  the  subroutine  TB  vhose 
listing  is  in  the  appendix.  The  subroutine  TR  calls  the  airfoil  subroutine 
after  determining  the  effective  Mach  number,  but  keeping  ell  other  parameters 
constant.  The  subroutine  TR  together  with  CORFAC  carries  out  the  steps 
listed  above. 


A  number  of  runs  have  been  made  of  these  programs  both  with  and  without  tip 
relief  to  determine  the  change  in  torque  resulting  from  tip  relief.  The 

o  +V  A  +V/>  AWO  A  ^  AAA  4  ^4  ..4  J4  Ai.4t.ai  J  -  J  _  J  1  .  .  i.  .  J  4  .. 

^  '-/AA  W44W  'AA.  A  A  UAO  UA  4UUVXWU  4.0  AilUl^a  VCU  4  U 

Figure  U.  As  expected,  the  drag  is  reduced  more  near  the  tip  and  aB  Mach 
number  increases.  The  largest  change  by  far  is  due  to  the  change  in  Mach 
number  rather  than  *tr  change  in  drag  coefficient  from  the  tip  relief 
theory.  A  further  ciucussion  of  these  results  is  postponed  for  a  few 
pages  so  that  comparison  can  be  made  with  a  simplified  strip  technique 
which  mny  be  used  in  energy  type  calculations. 


A  SIMPLIFIED  BLADE  ELEMENT  ANALYSIS  TO  CALCULATE  TIP  RELIEF 


By  making  a  simplifying  assumption  concerning  the  variation  of  the  drag 
coefficient  with  Mach  number  and  another  assumption  on  the  two  changes 
to  the  drag  coefficient  by  the  tip  relief  effect,  a  greatly  simplified 
equation  for  the  charge  in  rotor  torque  due  to  tip  relief  i t  found. 

The  result  of  these  assumptions  is  a  simplified  integral  for  torque 
due  to  tip  relief.  The  results  can  be  used  in  a  number  of  ways:  as  a 
general  indication  of  the  magnitude  of  the  torque  due  to  tip  relief; 
empirically  in  performance  calculations;  or  as  part,  of  an  energy  type 
computer  program,  since  computing  time  for  the  integration  routine  is 
very  short. 


In  a  previous  report1,  this  approach  vaB  described  for  the  hovering 
rotor.  Here  the  approach  is  extended  and  results  presented  for  a  rotor 
in  forward  flight. 


The  first  assumption  concerns  the  drag  coefficient  variation  vith  Mach 
number .  A  very  simplified  representation,  which  has  been  often  used, 
is  by  two  straight  lines  as 

CL  ■  CL  0  <  M  <  M'  (2U) 

D  D. 

inc 

C  -  ♦  a  (M  -  M’)  Ms  <  M  <  1.0  (2$) 

u  u.  c 

i  DC 


Here  M'  is  an  arbitrary  Mach  number  to  fit  the  approximation  and 
is  the  incompressible  drag  coefficient.  Since  will  drop  out  nC 

of  the  derivation,  its  value  and  variation  around  ?fie  rotor  disc  are 
ismaterial.  Typical  valueB  for  an  HACA  0012  airfoil  section  would  be 


M'  *  0.8 

a  «  0.57 
c 


(26) 


The  tip  relief  correction  to  the  drag  coefficient  is  made  up  of  two 
parts:  a  change  in  the  effective  Mach  number  resulting  in  a  change  in 
the  drag  coefficient,  and  a  change  of  the  drag  coefficient  itself. 

The  change  in  drag  coefficient  due  to  the  change  in  effective  Mach 
number  has  been  found  to  be  much  larger  than  the  change  of  the  drag 
coefficient  itself  in  the  subsonic  drag  rise  range.  Therefore,  the 
latter  ccrrecticr.  can  be  neglected.  By  neglecting  this  correction, 
there  is  no  tip  relief  effect  on  the  drag  coefficient  below  M* ,  since 
the  drag  coefficient  is  independent  of  Mach  number  in  this  range.  The 
equation  for  the  drag  coefficient  for  a  finite  wing,  in  the  range  M 
greater  than  M' ,  is  the  combination  of  the  two-dimensional  drag 
coefficient  with  the  tip  relief  correction 

CD  *  CD  +  ac  l(M  -  AM)  -  M'J  +  AC  (27) 

FW  inc  '  L' 


The  parameters  AM  and  AC^  are  the  tip  relief  corrections.  Neglecting 

AC  t>rm  rpflrrnn^i  ncr 

- -p  ^  * 


+  a  [M  -  M’  -  AM] 
DFW  Dinc  C 


(28) 


"d  cr.  +  cd 

Lnc  +  c  TR 


where 


Incompressible  two-dimensional  drag  coefficient 
of  the  section 


■  Drag  due  to  compressibility  effects  in  two-dimensional 

flnu  z  a  fM  _  M*  > 

"c . 
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*  Drag  due  to  the  tip  relief  effect  *  a  AM 


Using  equations  (19)  and  (20)  aod  Cjy^  as  defined  above  for  equation  (28), 


°D, 


TR 


a^AM  -  *r(l/0)(t/l*ir)(Ehn{6c/y)Io/TjM 


(29) 


Since  the  drag  coefficient  variation  with  Mach  number  is  taken  to  be  linear 
and  is  broken  into  its  components,  one  being  the  tip  relief  drag  coefficient, 
the  torque  coefficient  on  a  rotor  blade  may  be  similarly  broken  into  its 
components.  The  component  we  are  concerned  with  is  torque  due  to  tip  relief. 
With  the  usual  assumptions  concerning  small  angles,  no  radial  flow,  etc., 
the  torque  of  an  element  of  a  rotor  blade  due  to  drag  is 


dQ  ■  r  dD 


—  p(fir  +  V  Sin  <p)2  CD  2  c  r  dr 


(30) 


Since  the  equation  for  the  drag  coefficient  is  linear,  the  change  in 
torque  due  to  tip  relief  on  the  element  may  be  written  as 


dyTR  *  ~2  P  +  V  Sln  Cd_  2  c  r  dr 

TTv 

Transforming  to  coefficient  foi-n. 


(31) 


dC, 


_1  (fir  +  V  Sin  tl>)2  c  2  c  r  dr 


^TF  2  "  wF2p(nR)2E  ~TR 

with  Prandtl-Glauert  transformation  for  solidity 

o  =  b  2c  /itR  =  2b/"R  c ,  8  =  Bo. 
c  c  i  i 

and  the  definition  of  advance  ratio,  y 


dCQ  ■  (1/2)  (o  /be)  (M  Sin^+r)2  C  r  dr 
xR  1  UTB 


(32) 


( 33) 


(3M 


The  change  in  torque  per  blade  due  to  tip  relief  is  found  by  integrating 
the  elemental  torque  along  the  blade  span  and  averaging  around  the 
azimuth.  Mathematically  this  is  expressed  as 

Q  *  J  /R  dr  dy  (35) 

GO 


2* 


Thus 


CD 

iH 

o/b 


,  2n  i 

nr  /  /  f  ( u  sic  *  »  ~'2 


17  o  o 


+  r )  C_  r  dr  d(i 
TR 


(36) 
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How  Cjw.  io  given  by  equation  (29)  for  the  cue  where  only  the  6 H  effect 
of  tip^relief  is  used. 


^  1  2*  ]  _  r  t 

• —  ■  “  J  /  (u  Sin  i>  +  r)^  a  —7  X — 

o/b  4*  q  $  e  6*  4w 
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ec  i _  _  _ 

(Eh  { — }rtt>]  M  r  dr  d* 
n  y  T 


If  this  equation  is  restricted  to  hover  and  to  only  the  first  term  of 
the  Taylor  series,  it  vill  reduce  to  equation  (122)  of  Reference  1,  which 


is 

C 


OT& 


SrR 


I,  1 

l 


Mp 


(1  + 


Y-l 


i  - 


f, 


dr 


'  2&yXg ( 1-r ) 


♦  1 


(37a) 


7*he  integration  and  calculation  of  this  equation  are  done  by  the  use  of 
the  computer  programs  TKDRIVER,  TIPRLF  and  CORFAC.  (See  appendix.) 

The  programs  are  set  up  so  that  the  latter  two  subprograms  may  be  used 
with  an  energy  type  program,  while  TRDRIVER  will  drive  the  two 
subprograms  so  that  the  tip  relief  effect  can  be  evaluated  by  itself. 
The  program  TIPRW  carries  out  the  integration  using  a  rectangular 

fuxuiulu  us 

/  /  f(rlf*}  dFd*  -EE  f  (rj\  if  }  a?  a*  (38) 

0  0  i  j  i.J  i  J  i 

where 


AiJ/  *  2rr/n 


J  *  °«l  —  n* 


At  »  r  —  r 
i  i  i-1 


_  1  _ 
r .  -  (r,  +  r 


i-1 


)/2 


(39) 


and  the  values  of  rj  and  njj,  are  input  to  the  program.  The  tip  relief 
driver  has  three_sets  of  listed.  The  program  was  debugged  using  the 
first  set  of  10  r^'e.  Comparison  with  the  second  set  of  F^'s,  with  30 
values ,  indicated  a  big  difference,  but  comparison  of  the  third  set  of 
15  values  and  the  second  Bet  shoved  only  a  slight  difference.  Similarly, 
a  comparison  of  24  and  36  azimuthal  stations  did  not  show  any  significant 
differences. 


The  cost  of  computation  using  these  subroutines  is  minimal.  In  calcu¬ 
lating  the  plots  presented  here  on  a  DEC  PDP-10  computer,  some  serventy- 
five  ■‘"light  conditions  could  be  calculated  in  less  than  two  minutes  of 
CPU  time  at  a  cost  of  approximately  three  dollars  including  compilation. 
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In  the  previous  lections  of  this  chapter,  ve  have  described  a  theory  for 
the  quantitative  determination  of  the  tip  relief  and  tvo  calculative 
techniques.  In  this  section,  results  from  these  calculations  are 
shovn.  The  results  from  the  simplified  theory  are  presented  first,  as 
these  results  cover  a  wider  range  of  parameters.  The  number  of  calcu¬ 
lations  using  the  various  strip  methods  vas  limited  due  to  time  and 
cost}  th*  results  from  the  atrip  analyses  are  given  and  compared  to  the 
simplified  theory'. 

Ve  should  note  that  the  equation  derived  (see  equation  (37))  has  as 
the  result 


per  blade 


(o/b) 


which  is  identical  in  magnitude  to  the  quantity 


o 


uni.  the  first  quantity  indicates  tnai.  the  major  parameter  of  tip  relief 
is  solidity  per  blade  (o/b).  This  is  related  to  blade  aspect  ratio  by 

a.  jn  do) 

where  the  blade  aspect  ratio,  X,  is  defined  ss  the  rotor  radius  to 
chord  ratio.  Only  rectangular  planforms  are  considered  herein. 

Figures  5  and  6  show  the  variation  of  the  tip  relief  torque  with 
advancing  tip  Mach  number,  rotational  Mach  nvxaber,  and  advance  ratio 
aa  calculated  by  the  simplified  method.  The  two  figures  differ  in  blade 
aspect  ratio. 

These  figures  show  a  map  of  the  tip  relief  torque  variation  with  flight 
variables.  Lines  of  constant  rotational  Mach  number  (ftR/a)  would  be 
obtained  as  a  helicopter  gained  forward  speed  at  a  constant  rotational 
speed  (ft)  and  at  constant  temperature.  As  the  forward  speed  is  increased 
from  hover,  in  general  the  torque  due  to  tip  relief  increases  slowly. 

In  some  cases  of  low  (only  Blightly  above  M’ )  rotational  and  advancing 
tip  Macb  numbers  there  may  even  be  a  decrease.  The  slow  rise  or  decrease 
in  torque  occurs  because  as  the  Mach  number  decreases  on  the  retreating 
side,  the  tip  relief  effect  is  also  reduced,  vhile  there  is  the  increase 
in  Mach  number  and  tip  relief  torque  on  the  advancing  side.  Their 
balance  determines  the  magnitude  of  the  decrease  and  the  rate  of  rise 
of  tip  relief  torque  at  these  low  forward  speeds.  After  the  Mach  number 
on  the  retreating  Bide  has  decreased  to  a  negligible  value,  the 
IV  muuw  1 — t  vjyc  ci  A  cv  i#  uct.  lou  CD  Mfl  i  &OU  t-  flBCM  • 


21 


At  very  moderate  advance  ratios,  above  about  0.10  to  0.1^>,  the  curves 
coalesce  and  become  independent  of  either  rotational  Mach  number  or 
advance  ratio.  A  line  indicating  the  lover  envelope  of  these  points 
has  been  drawn,  which  is  vsed  in  the  following  comparisons ,  and  is 
labeled  "high-speed  envelope". 

Thus,  as  far  as  tip  relief  is  concerned,  at  rotational  tip  Mach  numbers 
below  M* ,  the  value  of  the  tip  relief  torque  depends  mainly  on  advancing 
tip  Mach  number  and  only  slightly  (for  practical  purposes,  negligibly) 
on  advance  ratio  or  SlR/a. 

Figure  7  indicates  how  tip  relief  varies  with  various  parameters  and 
compares  calculative  techniques  and  some  limited  test  data.  This  is 
discussed  in  detail  in  the  following  paragraphs. 

The  extremes  of  these  two  sets  of  curves  from  Figures  5  and  6,  the 
hover  xnd  the  high-speed  envelope,  are  replotted  in  Figure  7(a)  to 
■how  the  effect  of  solidity  per  blade  on  the  tip  relief  torque.  The 
tip  relief  is  calculated  for  a  semi-infinite  wing  tip  relief  model, 
it  does  depend,  to  am  important  degree,  on  aspect  ratio. 

The  key  parameter  for  tip  relief  is  distance  from  the  tip  in  terms 
of  chord  and  local  Mach  number.  The  tip  relief  effect  diminishes  rapidly 
with  distance  from  tip.  The  rotor  blade  with  the  greater  chord  (i.e., 
lower  aspect_ratio,  a  higher  solidity  per  blade),  will  have  an  element 
at  the  same  r,  and  thus  the  seme  Mach  number,  closer  to  the  tip  ip 
terms  of  chord  than  a  blade  with  a  shorter  chord.  This  element  will 
therefore  have  a  larger  tip  relief  correction.  Thus  the  blade  with 
the  greater  chord  will  have  greatest  tip  relief  and  the  larger  allevia¬ 
tion  of  compressible  power  loss,  other  things  bting  equal. 

Figures  7(a),  (b)  and  (c)  present  comparisons  between  the  simplified 
calculation,  atrip  analysis  and  some  flight  test  data.  These  comparisons 
must  be  analysed  with  care,  as  all  of  the  methods  have  certain  problems 
associated  with  them.  The  simplified  calculations  ta'ie  too  few  points 
near  tip  Mach  numbers  of  0.60  and  diverge  greatly  near  Mach  one  due  to 
the  Prandtl-Glauert  effect.  The  two-dimensional  airfoil  data  at  high 
subsonic  Mach  numbers  may  not  be  as  valid  as  one  would  like. 

In  hover  (Figure  7(b)),  the  simplified  theory  gives  a  much  higher  value 
of  tip  relief  than  does  strip  analysis.  On  the  other  hand,  at  a  high 
advance  ratio,  the  agreement,  as  shown  in  Figure  7(c),  could  not  be 
better.  Ho  explanation  for  this  discrepancy  and  this  agreement  has 
been  found.  But  note  thst  computation  for  both  cases  vas  made  only  in 
the  range  of  0.82  to  0.09  in  advancing  tip  Mach  numbers. 

Presented  in  Figure  7(d)  is  some  comparison  with  test  data.  The  method 
of  determining  the  experimental  values  of  torque  due  to  tip  relief 
requires  some  detailed  explanation.  Reference  3^  (specifically  in 
Figure  9)  presents  both  test  data  and  theoretical  calculations  for  the 
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per  blade 


Advancing  Tip  Mach  Busbar 


Figure  5.  Change  in  Torque  Coefficient  I>ue  to  Tip  Relief  (o/b  0.0159). 
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Test  and  Predicted 
Power  Req-iired .  34 


the  test  data  increases  with  ndvu: 
the  theory  v-uld  predict. 


variation  of  the  power  coefficient 
with  advancing  tip  Mach  number  for 
the  l  i-lF  helicopter  at  p  *  0.25. 

The  test  data  is  presented  as  a 
curve  (as  opposed  to  data  points) 
and  includes  instrument  corrections, 
cross  plottiixg,  fairing,  etc. 
Performance  prediction  ia  based  on 
Tanner18  and  includes  the  volidity 
coi ruction.  The  Pull  difference 
between  i  he  power  required  ^est  data 
of  Reference  0  and  the  'alculated 
value  from  the  Tanner  charts26  is 
attributed  to  the  tip  relief  effect, 
A  typical  curve3*1  and  the  power 
attributed  to  tip  relief  are 
sketched  in  Figure  8.  binee  ir 
magnitude  Cp  ~  Cq,  the  value  of  Cq^ 
determined  in  this  catuie."  is 
divided  by  a  and  replotted  in 
Figure  ?(d),  it  is  surprising  that 
icing  tin  Much  number  more  rapidly  tiian 


It  shouDd  be  noted  in  regard  to  the  accuracy  required  tr.at  the  calculation 
is  of  the  change  in  torque  coefficient  due  to  tip  relief.  The  order  of 
magnitude  of  Cq  is  about  25  x  10~8,  while  at  the  same  time  the  charge  in 
torque  coefficient  due  to  tip  relief  is  about  2  x  10~5.  Thus,  there 
is  a  raptor  of  1C  between  these  two  numbers,  which  are  subtracted  from 
each  other.  Three-place  accuracy  or  &Cp~,  requires  four-place  accuracy 
of  Cq. 

The  results  of  the  theory  overall  appeal  to  be  useful  in  performance 
calculations.  Some  more  detailed  calculations  are  needed  to  resolve 
some  anomalies,  especially  calculations  using  strip  theories.  Also, 
further  comparisons  with  test  data  would  be  useful. 

rr<rp  uri  i  re  COMBI  NET’  WITH  COMPRESS  I  BTI.I  TV 


It  should  be  emphasised  that  the  calculation  of  the  tip  relief  as  pre¬ 
sented  here  is  not  a  calculation  of  the  effect  of  compressibility;  rather, 
it  is  a  correction  to  the  calculation  of  compressibility  effects,  au 
presently  used  for  rotor  performance  calculations  using  two-dimensional 
airfoil  data.  In  strip  analysis,  tip  relief  corrects  the  two-dimensional 
compressible  flow  airfoil  coefficients.  In  energy  type  calculations, 
core  must  be  taker,  to  be  sure  that  tip  relief  is  not  already  included  in 
the  method  used  to  calculate  coi-pr essibil ity  power  losses.  For  example, 
empirical  curves  fitted  to  tcr-t  data  will  by  definition  include  the  tip 
relief  effect. 

To  indicate  the  importance  cf  this  point ,  the  calculation  cf  power  loos 
due  to  compressibility  1  rom  Reference  13 2  has  bf'-.,  combined  with  the 


/simplified  tip  relief  ct4  u  n’t.'--':',  as  presented  here.  Gtssov  and  Crinr*- 
calculated  the  compress; Vie  pro.-.-  loss  using  a  strip  analysis-1 7  with 
compressible  airfoil  dat.\.  The  compressible  power  loss  is  defined  as  the 
difference  betwee-n  the  power  at  a  lev  and  a  high  tip  speed.  Figure  9 
indicates  the  results  of  Gessov  and  Crim32  later  calculations  presented 
by  Dingeidein,45  flight  test  results,33  arid  the  combined  compressibility 
and  tip  relief  power  loss. 

The  key  result  of  calculating  the  compressibility  power  loss  of  Gessow 
and  Crim3^  for  a  series  of  edvance  ratios  (from  0.2  to  0.5).  twists 
(0  *  and  — l6  ),  and  solidities  (0.025  to  0.08)  was  that  the  most 
important  parameter  is  the  differercr.  between  the  advancing  tip  Mach 
number  and  the  drag  divergence  Mach  number  (tip  ) .  They  drew  a 

single,  best-f 1 t  curve  through  the  celeulated  point  on  a  Cp ?/o  vs. 

Mi .90  *  Mdd  Plot.  This  is  the  curve  labeled  in  Figure  9  as  being  from 
Reference  32.  We  should  note  that  the  scatter  is  of  the  order  of 
±0.001  ir  Cpc/o,  which  is  one  division  on  +he  ordinate  scale.  Reference 
t*5  found  an  even  lsrger  spread,  as  indicate!  by  the  two  limit  curves  in 
Figure  9. 

The  high-speed  envelope  value  of  Cq^/.i  is  subtracted  from  that  of  Gessow 
and  Crim  2  to  obtain  the  tip  relief  corrected  compressibility  power 
loss  curves.  Mean  values  of  Cq  /a  from  the  two  solidities  (0.0159 
and  0.0232  of  Figures  5  and  6)  are  used,  as  the  other  curves 
represent  a  wide  range  of  solidities.  For  the  first  curve  we  used 

M 1  =  M  -  «  =  n  Wo  ( rs«v-*r*.  -i  \  n..+  u.  ^  ^ 

•*uu  ^  ^  A-*aa  .'.u  i\a  ttt  tructr  2d  ta  fcfcl'O  fcurigi.e  Ol 

attack  is  about  C.74;  therefore,  Cq^/o  was  also  calculated  with 

M'  =  0.74  and  using  ac  =  0.492  to  represent  the  rirfoil  (Curve  2). 

Hie  tip  relief  corrected  compressibility  power  loss  becomes  almost  a 
straight  line  in  both  cases  when  plotted  against  Mp  cL)  -  Kdd  (Figure  9). 
The  Prandtl-Glauert  law  increases  the  tip  relief  at’the  same  rate  that 
compressibility  increases  the  power  loss  ,  When  compared  to  the  AFFTC 
flight  test,  results,33  the  prediction  is  still  greater  than  the  test 
values  which  would  be  expected  co  include  tip  relief.  For  a  more  com¬ 
plete  comparison,  a  better  representation  of  the  drag  coefficiem  than 
equation  (25)  is  needed. 
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Figure  9*  Comparison  of  Compressibility  Effects, 

TIP  RELIEF  CALCULATION  AS  PF0P3S iP  BY  PROU-"’. 

Theory  -  Description,  Discussion  and  Critique 

In  the  October  1971  issue  cf  the  AHS  Journal.  P^outy46  described  a 
simple,  empirical  means  of  calculating  the  magnitude  of  the  tip  relief 
effect.  The  approach  states  that  the  apparent  three-dimensional  drag 
divergence  Macli  number  on  a  rotor  differs  from  the  tve -dimensional 
measured  drag  divergence  Mach  number  of  a  "correction"  tfhich  vas 
previously  developed  by  Sophnr.5  A  farther  empirical  multiplication 
factor  for  the  ferward-fligh-  case  accounts  for  the  fact  that  "until  a 
significant  portion  of  the  blade  dvells  in  a  region  of  dr&%  divergence 
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for  a  significant  time  the  effects  are  not  apparent  in  measurements 
of  power." 

This  is  similar  to  the  simplified  method  presented  here,  in  that  a 
correction  to  the  Mach  msaber  is  made.  The  nature  of  the  blade  element 
integration  illl  account  for  the  second  factor. 

An  examination  of  Trouty's  approach  to  tip  relief  appears  to  indicate 
that  it  is  based  on  a  theory  which  is  not  applicable  to  the  problem. 

In  the  following  paragraphs,  the  Mach  number  correction  as  proposed 
by  Prouty  will  be  described  along  with  its  theoretical  basis  and  the 
reasons  for  doubting  its  theoretical  validity. 


Prouty  applies  the  correction  to  the  drag  divergence  Mach  numbe",  but 
he  indicates  that  the  correction  can  be  applied  to  other  than  this 
particular  value  of  the  Mach  number.  The  Prouty  correction  factor  is 


(hi) 


where  Cp^  is  an  average  incompressible-flow  airfoil  pressure 
coefficient.  For  the  basis  of  using  this  correction  factor,  Prouty 
simply  states  that  it  was  derived  in  Sopher's  work,5  who  in  turn  refers 
to  KUchemann  and  Weber.47 


A  relationship  between  the  theoretical  potential  flow  equations  for 
compressible  flow  and  those  for  incompressible  flow  is  the  Prondtl- 
Glauert  transformation 


P2  =  1  ~  (h2) 

A  further  refinement  for  a  transformation  between  theoretical  com¬ 
pressible  flow  and  theoretical  incompressible  flow  was  developed  by 
KUchemann  and  Weber47  and  is  given  as 

e2  =  1  -  M*  (l  -  cp  )  (h3) 

"  inc 

In  the  original  paper,47  it  was  suggested  that  C-pi  should  be  a 
suitable  chosen  mean  value  cl  the  pressure  coeffic?ent  over  the 
suction  region  of  the  airfoil  in  incompressible  flow.  Later  similar 
extensions  of  the  Frand  -,1-Glauert  transformation  all  use  the  incom¬ 
pressible  preesure  coefficient  but  suggest  that  the  local  rather  than 
a  mean  value  be  used  (for  example.  Lock,  et  al48). 

The  extension  of  the  KUchemann -Weber  rule  to  rotors  was  made  b>  Sopber5 
by  use  cf  f.  to-called  "equivalent  linear  theory  tip  Mach  number."  If  Cp, 
is  e.  constant,  then  the  quantity  M^l  -  ^P^nc  can  be  considered  as  c 
a-  "equivalent"  Macb  number.  In  the  case  of  a  rotor,  M»  is  replaced 
V/  illi/t.,  and  thus  (fiR/a)/T’ -  ^rlQC  would  represent  an  "cquivaleul 
lirtar  theory  tip  Me.ch  number.  ' 
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Equation  (*»3)  was  derived  for  compressible,  two-  and  three-dimensional 
flew  over  airfoils  and  wings,  and  is  the  transformation  between  theoretical 
compressible  flows  and  an  equivalent  incompressible  flow.  It  is  net 
related  in  any  way  tfcnt  ve  can  determine  to  any  relation  between  tvo- 
and  three-dimensional  flovB. 

In  conclusion,  the  applicability  of  the  theoretical  basis  for  Prout.y's 
approach  to  tip  relief  is  doubtful. 


DISCUSSION  OF  SIMILARITY  LAWS  FOR  ROTORS 


The  concept  of  similarity  lava  is  important  for  the  application  of  fluid 
mechanics  in  design  practice.  The  governing  differential  equation  of 
fluid  motion  lends  itself  to  a  solution  only  in  rare,  simplified 
instances.  But  an  examination  of  these  equations  and  the  boundary  con¬ 
ditions  can  lead  to  similarity  lavs  which  can  then  be  used  to  assess 
the  effects  o i  varying  various  characteristics  of  the  flow. 


Similarity  lavs  group  the  parameters  (including  nondimens ional  parameters) 
of  the  flow  into  functional  groups,  while  dimensional  analysis  gives 
only  the  nondimens ional  parameters  such  as  thickness  ratio,  Mach  number, 
etc.  The  best  knovra  of  the  similarity  lavs  is  the  Prandtl-Glauert  lav. 
Because  the  Prandtl-Glauert  lav  is  derived  only  on  the  basiB  of  a 
linearized,  homogeneous  equation,  it  may  be  stated  in  a  number  of  ways, 
for  example. 


C 


P 


r 


(kk) 


for  a  two-dimensional  airfoil  section.  In  this  form,  the  Prandtl-Glauert 
lav  states  that  the  pressure  distribution  vill  be  the  same  for  a  family 
of  airfoils  of  different  thicknesses  at  different  Mach  numbers  if  the 
quantity  t / /l-Mz  is  kept  constant.  Hie  pressure  distribution,  of  course, 
when  properly  integrated,  gives  the  drag,  lift  and  moment;  therefore, 
the  functional  variable  for  drag,  lift  and  moment  is  very  similar  to 
the  form  given  in  equation  (kh). 

The  speed  range  of  interest  for  high-speed  helicopter  rotors  is  the  high 
subsonic  and  transonic  region.  Similarity  lavs  for  transonic  flows 
were  developed  by  a  number  of  investigators;  a  complete  description  and 
bibliography  is  given  by  Spreitei  .  The  transonic  similarity  rule  for 
drag  is 

CD  =  CD  (M,X }  (U5) 


where 


M 


T5/' 

.  _ _ 

[m2(y+i)t]2/3 


X  =  [M2(y+1)t J1/3  X 


The  quantities  given  by  C^,  M  and  X  are  called  reduced  drag  coefficient, 
reduced  Mach  number,  and  reduced  aspect  ratio.  Because  of  the  functional 
nature  of  the  similarity  rule,  t  is  relation  may  be  given  in  somewhat 
different  forms;  The  similarity  rule  simply  states  that  a  plot  of 
against  either  M  or  X,  while  holding  the  ether  constant,  vill  result 
in  a  single  curve.  The  similarity  rule  by  itself  does  not  give  the 
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equation  for  this  curve.  It  nay  be  possible  to  empirically  cross-plot 
airfoil  dc la  and,  for  certain  ranges ,  obtain  an  equation  relating  the 
reduced  variables. 

The  tip  relief  theory  for  either  fixed  wings2  or  rotary  wings1  is  a 
similarity  law,  in  that  they  give  the  relationship  of  the  drag 
variation  between  two-  and  three-dimensional  ad r foils.  The  theory 
as  derived  is  not  in  a  functional  fort,  hut  in  application  the  tip 
relief  theory  does  reduce  the  number  of  curves  required  to  present  the 
Information.  For  example,  to  present  the  drag  coefficient  of  wings  of 
different  aspect  ratios,  a  separute  curve  would  be  required  for  each 
aspect  ratio.  By  applying  the  correction  factor  from  the  theory,  this 
same  6et.  of  data  can  be  presented  as  a  single  curve.2  This  differs 
from  the  usual  application  of  similarity  rules  where  the  data  is  pre¬ 
sented  with  a  new  set  of  coordinates. 


The  possibility  of  comparing  the  results  of  the  similarity  rule  given 
by  the  tip  relief  theory  and  the  conventional  transonic  similarity  rule 
has  been  investigated  without  any  success  to  date.  A  notable  feature 
of  the  transonic  similarity  rule  relating  drag  coefficient  to  aspect 
ratio  i6  that  thickness  is  raised  to  powers  other  than  the  first  in 
the  functional  variable,  as  is  shown  in  the  example  for  drag  in  equation 
(45).  On  the  other  hand,  the  tip  relief  theory,  based  as  it  is  on 
linearized  airfoil  theory  and  the  linearized,  elliptic  equation  of 

OllVte  /“»  ;  «-<*.♦  f  rfinffn  r>.A  4-  ~  -v _  _  _  -f  V  *1  - 

gives  a  variation  of  thickness  ratio  to  the  first  power. 


lUlit)  f 


Further  work  in  applying  transonic  similarity  lavB  to  helicopter  rotor 
aerodynamics  is  needed.  Derivation  of  similarity  rules  for  rotorB, 
such  as  those  available  for  fixed  wingB ,  would  be  very  useful,  but 
it  appears  to  be  more  in  the  vay  of  a  long-term  goal.  More  immediate 
applications  would  include  tests  of  the  validity  of  the  airfoil  data 
nov  in  use  and  their  extension  to  higher  speeds,  and  the  use  of 
similarity  rules  to  describe  airfoil  data  in  energy  type  calculations . 
Thia  latter  application  would  extend  the  present  use  of  tvo-dimeD3ional 
data  in  rotor  aerodynamics  by  use  of  the  transonic  similarity  ruleB. 
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COMMHITS  OH  EXPERIMENTAL  VISIPICATIOB  OP 
COMPRESSIBILITf  EFFECTS  AHD  TIP  BEUff 


There  is  solid  theoretical  evidence  of  the  existence  of  a  compressibility 
tip  relief  effect,  on  helicopter  rotors.  There  are  similar  experimental 
trends  to  that  effect.  Assuming  the  theoretical  model  to  be  sound, 
it  vould  be  desirable  to  design  the  definitive  experiment  to  ’'calibrate" 
tip  relief.  This  presents  great  difficulties,  vhich  vill  be  discussed 
here. 

Tests  for  the  overall  tip  relief  effect  for  fixed  vingo  in  uniform 
flov  were  described  in  Reference  50,  where  very  good  correlation  was 
found.  In  the  extension1  to  rotary  wings,  two  steps  are  used:  First 
the  variation  of  the  tip  relief  correction  along  the  span  is  calculated 
for  uniform  flow.  The  second  step  combines  this  variation  with  blade 
element  theory;  namely,  the  tip  relief  correction  of  a  particular 
element  in  a  rotor  is  the  same  as  the  tip  relief  of  the  element  in  a 
fixed  wing  that  is  located  at  the  same  distance  from  the  tip  in  a  flow 
having  the  free-stream  velocity  of  the  blade  element.  The  Mach  number, 
Reynolds  mstber,  angle  of  attack,  etc.,  are  also  the  same  for  the  tvo 
elements.  Thus  there  are  tvo  hypotheses  that  need  experimental 
verification. 

Testa  of  the  spanwise  variation  of  the  tip  relief  effect  could  be  per¬ 
formed  in  a  very  large  transonic  wind  tunnel,  since  the  effect  is  mainly 
in  the  first  chord  length  near  the  tip,  one  would  need  the  detailed 
drag  distribution  In  this  region.  The  need  for  a  large  chord  so  as  to 
be  compatible  with  the  instrumentation,  followed  by  the  need  for  a 
reasonable  aspect  ratio  and  a  large  throat  to  avoid  interference,  would 
require  an  enormous  wind  tunnel.  As  a  beginning  in  testing  for  the  tip 
relief  effect,  local  pleasure  distributions  vould  be  useful.  Systematic 
measurements  of  the  pressure  distribution  on  wings  of  various  aspect 
ratios  at  high  Mach  numbers  are  needed.  At  incompressible  speeds, 
some  data  of  this  nature  is  available,51  but  the  change  due  to  tip 
relief  ia  less  than  the  accuracy  of  the  experiment,  and  comparison  is 
not  possible.  Compressible  airfoil  data  appear  to  be  mostly  measured 
by  balances.52*53 

Tests  of  the  application  to  rotors,  by  systematic  variation  of  parameters 
which  effect  the  tip  relief  vitbout  changing  other  aspects  of  rotor 
performance,  would  be  more  intellectually  satisfying. 

As  stated  several  times  in  this  report,  tip  relief  is  correction  to  the 
preoent  calculation  of  compressibility  effects.  Assuming  that  one  can 
differentiate  in  a  given  set  of  test  data  between  incompressible  and 
compressible  pover,  one  must  further  be  able  to  fully  account  for 
ccangea  in  two-dimensional  compressibility  effects  before  one  can  calcu¬ 
late  the  three-dimensional  effects,  i.e.,  tip  relief.  Today,  for 
example,  reliable  tvo-dimensional  test  results  for  the  0012  airfoil 
at  Mc,Cu  mU.x.o  i'j  uuc  cue  uuL  y.si,  avuliauie .  It  can  be  said  unau 
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detailed  information  on  helicopter  rotor  compressibility  effects  is 
sketchy  at  test.  The  determination  of  tip  relief  therefore  requires 
trustworthy  calculation  of  the  rotor  performance  baaed  on  blade  element 
analysis . 

Since  the  tip  relief  compressibility  effects  are  localized  near  the  tip, 
it  is  tempting  to  try  to  test  a  blade  tip  element  in  a  stationary 
position,  vith  a  spanwise- varying  simulated  flow,  at  transonic  speeds. 

We  were  not  able  to  cone  up  with  a  reasonable  wind-tunnel  configuration 
to  accomplish  this  objective,  and  we  conclude  that  any  tip  relief  tests 
must  be  made  with  a  rotating  blade  and  a  constant  free-stresm  flow,  or 
in  hover. 

In  this  report,  limited  comparison  to  flight  test  data  is  made 
(Figure  7(d)).  There  is  an  increasing  amount  of  flight  test  data  at 
compressible  tip  speeds  becoming  available,  and  Blmilar  comparisons 
of  these  data  are  recommended.  Time  did  not  allow  these  comparisons  to 
be  performed,  because  no  time  was  allotted  to  calculate  the  performance 
of  the  available  flight  te6t  cases  using  present  performance  prediction 
methods.  The  experimental  value  of  tip  relief  is  taken  as  the 
difference  between  the  calculated  prediction  and  the  measured  power 
required.  This  leaves  possibilities  for  errors  in  a  number  of  ways, 
including  the  estimate  of  fuselage  drag,  and  the  effects  of  variable 
induced  velocity,  blade  flexibility  and  unsteady  effects.  Thus  in  the 
flight  test  data,  differences  due  to  other  causes  may  be  thought  of  as 
errors  in  the  tip  relief  theory. 

In  principle,  one  can  verify  tip  relief  in  a  manner  similar  to  that  used 
by  Anderson2 »50  in  the  fixed-wing  case.  This  requires  the  use  of  three 
or  more  physically  different  rotors.  Each  rotor  is  identified  by  its 
blade  aspect  ratio,  or  o/b.  Preferably,  o/b  should  vary  as  much  as 
possible  for  the  three  rotors,  say,  at  least  by  a  factor  of  2  tc  1.  One 
then  measures  (and  calculates  from  strip  analysis )  C  /a  against  90 
for  various  Mach  numbers,  say,  between  0.7  and  0.9,  at  a  constant  &r/o. 
Assuming  now  that  one  has  three  rotors,  identical  except  for  the 
parameter  o/b,  each  has  a  different  Cq/o  against  M-^  curve,  for  a 
given  C'  .  However,  these  curves  should  superimpose ’if  the  tip  relief 
correction  of  this  report  is  applied.  This  method  can  be  used  either 
for  wind-tunnel  teBts  or  for  flight  tests. 


C0HCLU3IQHS  AMD  RKXMffiNDATIOSS 


The  methods  used  to  calculate  the  effects  of  compressible  flow  on  the  per¬ 
formance  of  helicopter  rotors  have  been  reviewed.  The  application  of  a 
theory  for  tip  relief  based  od  the  complementary  wing  concept  has  been 
described,  and  results  from  various  applications  have  been  given. 

The  results  of  the  calculation  and  tip  relief  by  the  complementary  wing 
theory  Indicate  that  the  technique  can  and  should  be  included  in  blade 
element  type  rotor  programs.  This  is  done  by  modifying  the  two-dimensional 
airfoil  data  now  being  used  by  the  technique  presented  here  to  obtain 
three-dimensional  airfoil  data,  which  better  represents  the  actual  three- 
dimensional  flow.  Simplification  of  the  strip  analysis  and  the  tip  relief 
theory  allows  the  calculation  of  the  power  associated  with  tip  relief  using 
a  very  quick  computer  program.  These  resultB  match  reasonably  well  with 
both  tests  and  the  more  complex  strip  programs \  thus  this  simplified  tech¬ 
nique  should  also  be  included  in  energy  type  calculations  when  combined 
with  proper  representation  of  the  compressibility  power  required. 

More  work  Is  needed  in  this  area  of  properly  predicting  the  performance  of 
rotors  operating  in  the  compressible  flow  regime.  There  iB  available  a 
large  body  of  test  data,  including  the  work  of  Carpenter  at  NACA-Langley , 
the  Edwards  flight  teste,  and  the  NASA  Ames  wind  tunnel  tests,  and  a 
systematic  comparison  of  the  various  prediction  techniques  reviewed  here 
with  this  rang"'  of  test  data  is  needed.  By  making  such  a  systematic  com¬ 
parison,  the  best  of  the  prediction  techniques  can  be  chosen  and  developed 
for  use  in  strip  and  energy  calculations. 
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APPENDIX  -  COMPUTE  PROGRAMS 


The  four  computer  programs  d.'scusstd  in  this  report  are  listed  in  tne 
following  pages.  These  programs  are  short  and  very  straightforward.  Thus 
a  detailed,  description  will  not  be  given. 


CORF.aC 


Th  -S  program  calculates  the  value  of  the  right-hand  side  cf  equation  ( 19 )  • 


Call  Statement  Variables 


D'JOU 

- 

the  correction  factor  =  value  of  equation  (19) 

TAU 

- 

thickness  ratio 

BETA 

- 

Prandtl-Glauert  variable 

MACH 

- 

Mach  number 

RBA 

- 

radial  station  as  fraction  of  radius 

RAD 

- 

radius  in  feet 

CHD 

- 

chord  divided  by  diameter 

Switches 

1FR 

if  equal  to  zero,  does  not  print;  it  greater  than 

zero,  prints 

IFW 

CORFAC  variables 

if  equal  to  one,  uses  finite  wing;  if  equal  to  zero,  uses 

NTFRM 

* 

semi- infinite  vir.g* 

if  equal  to  three,  uses  three  terms  of  the  Taylor 

series** 

TP 

This  subprogram  calculates  the  dreg  coefficient  CD  using  equations  (20) 
through  (23),  the  subroutine  CORFAC,  Rnd  an  airfoil  subroutine  ERFOYL.  Most 
of  the  variables  pertain  to  ERFOYL  and  thus  will  not  be  described. 

Switches 

IMD  -  if  equal  to  one,  only  Mach  correction  is  calculated 

IPRIirr  -  if  not  c<4uui  to  zero,  then  results  at  this  point  are  printed. 

TIPRLF 

The  subroutine  TIPRLF,  using  CORFAC,  generates  a  value  of  the  change  in 
torque  coefficient  per  blade  and  the  change  in  torque  coefficient  divided 
by  solidity  due  to  tip  relief.  CORFAC  provides  the  tip  relief  correction 
factor  while  TIPRLF  handles  the  integration  process.  The  two  subroutines 
are  designed  for  final  use  with  energy  type  programs.  The  integration  uses 
the  value  of  the  correction  factor  at  the  mean  of  the  interval,  namely  at 


ri  +  rj-l 


*  As  listed,  IFW  ■  0 
»'•*  As  listed,  NTERM  =  3 

Ul 


Program-Sub program  Relationship 


Identification: 

Call  TIPRLF  'XNDA,  OMC-R,  COSAL,  VFPS ,  Cl,  XML1 ,  NPSI,  AC,  DCQTR ,  RAJ),  CIID, 
TAU,  SIGPB ,  BN,  D'JQTRS ,  NSTA,  NR  ST ,  I  PR) 

COMMON 


The  values  of  the  radial  station  to  be  used  are  in  the  COMMON  statement. 
Description: 

The  integration  of  a  simplified  blade  element  analysis  for  the  tip  relief 
correction  is  done  using  the  mean  value  within  the  interval.  Cir.ce  the  pro¬ 
gram  is  straightforward  and  not  lengthy,  only  a  few  exceptions  are 
described.  The  value  of  the  integrand  is  calculated  and  summed,  along  the 
radius  before  advancing  to  a  new  azimuth.  If  the  advance  ratio  is  less  than 
0,001,  only  one  azimuthal  station  is  calculated,  and  hover  is  assumed. 


Warning  Statements: 


b. 


V  v ~  ’  n  one 

“(1,90) . 


ou  VV1U0< 


ADVANCING  TIP  MACH  NUMBER  EXCEEDS  ONf 


*-  V  . -  ~  3. 

aw  u  t  mo  ±Ji.  j.  a  wt  u  • 


:i,9o) 


,=>x.yz 


If  the  local  Mach  number  at  the  first  radial  station  exce.0" 
the  input  value  of  the  drag  divergence  Mach  number,  the  it nr 
is  printed: 


SPAR  IS  ABOVE  DRAG  DIVERGENCE  MACH  NUMBER 
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MNEMONICS  LIST 
Subroutine  TIPRLF 


Program 

Mathematie.-d. 

Explanation 

Symbol 

Model  Rare 

AC 

a 

c 

Slope  of  compressible  drag  rise 

BN 

b 

Number  of  blades 

BETA 

6  =  /Tm* 

Prandtl-Glaucrt  factor 

Cl 

a 

Speed  of  sound 

CUB 

c/2R 

Chord  diameter  ratio 

COSAL 

con  a 

Cosine  rotor  angle  of  attack 

DCQTm 

C 

QTR 

Change  in  torque  coefficient  for  BN  blades 

UCQTRS 

C  ^,.,/c 

Change  in  torque  coefficient  divided  by 
solidity  for  the  BN  blades 

DELPSI 

A\p 

Step  in  azimuthal  angle 

DUOU 

AU/U  =  G 

«  X 

Correction  factor 

GAM 

r  =  1  +  m2 

Stagnation  compressibility  factor 

IPR 

Switch  for  detail  printing 

M190 

M1 ,90 

Advancing  tip  Mach  number, Note  1 

ML 

M 

Local  Mach  number, Note  1 

1CLT 

M 

‘  l,* 

Tip  Mach  number  at  azimuth, Note  1 

NF5I 

360/ Aijj 

Number  of  azimuthal  stations 

NRST 

Largest  number  of  radial  stations  used  for 
the  integration  at  any  azimuth 

NST 

Total  number  of  stations  used  for  the 
integration 

OMG.K 

hi: 

Rotational  tip  speed  in  feet  per  second 

PI 

=  3.1^16 

PS  I 

Azimuth  angle,  radians 

R 

r 

Raaial  station  as  fraction  of  radius 

Note  1. 

These  valuables  ore 

declared  as  REAL. 

43 


Explanation 


Program  Mathematical 


RAD 

R 

Rotor  radius 

SIGPB 

a  /b 

Solidity  per  blade 

SIHPSI 

Sin  ^ 

SUM 

I 

Summation  of  integrand 

TAU 

X 

Airfoil  thickneBa  ratio 

UL 

*=  p  Sin  $  +  r 

Local  normal  velocity 

VALUE 

Integrand 

VFPS 

V 

Forward  apeed 

XMDA 

Drag  divergence  Mach  number 

XMU 

P 

Advance  ratio 

TIP  RELIEF  CALCULATION  DRIVER 

This  la  a  simple  driver  for  TIPRLF  written  so  that  a  particular  rotor  con¬ 
figuration  is  inputted  and  tae  tip  relief  torque  coefficient  i*  calculated 
at  a  aeries  of  advance  ratios. 

Input 

The  input  data  is  put  on  three  unformatted  cards: 

Card  1 


RAD,  BN,  0(4®,  CHD 

Card  C 


TAU  ,  XMDA,  Cl,  AC,  CORAL 
Card  3 


LCASE,  IPR.  NPSI ,  LXMU,  UXMU,  DXMLJ,  WRSET 


The  variables  names  on 
T1PRLF. 

the  first  two  cards  are  identical  to 

thoae 

in 

Variable  Name 

Explanation 

Card 

Position 

Rotor  Characteristics : 

AC 

slope  of  compressible  drag  rise 

2 

U 

BN 

number  of  blades 

1 

2 

CHD 

chord  diameter  ratio 

1 

U 

COSAL 

cosine  of  rotor  angle  of  attack 

2 

5 

Cl 

speed  of  sound 

2 

3 

CWCR 

rotational  tip  speed 

1 

3 

RAD 

radius 

1 

1 

TAU 

airfoil  thickness  ratio 

2 

1 

XMDA 

drag  divergence  Mach  number 

2 

2 

Program  Switches 

IJMU  (Declared  Real) 

UXMU 

lower  (  used  in  determining 

U?*>er  1  u’s  in  the  calculations 

3 

3 

U 

5 

6 

DXHJ 

etep  V. 

3 

LCASE 

case  number,  negative  for  last  case 

3 

1 

IPR 

detailed  printing  if  positive 

3 

2 

HP8I 

number  of  azimuthal  stations 

3 

3 

IRSET 

number  (label)  of  set  of  radial 

3 

7 

stations 


t 


i 


Restrictions 


Up  to  three  sets  of  radial  stations  may  be  used,  each  with  up  to  30 
stations. 


COMOfi 

The  values  of  the  set  of  radial  stations  to  be  used  (determined  by  NRSET) 
are  transferred  by  the  CO>MON. 

Block  Data: 

The  values  of  the  three  sets  of  radial  stations  are  entered  through  block 
data  statements. 


Output 


All  output  variables,  except  with  detail  print,  are  labeled.  First  all 
rotor  characteristic  variables  are  printed.  Then  a  tabulation  of  p,  v. 


C<JTR/° 


is  given.  The  radial  stations  used  are  printed  below 


Note : 

After  a  tip  Mach  number  of  one  is  reached,  only  two  more  values  of  the 
advance  ratios  are  used  for  calculations  before  going  to  the  next  set  of 
inputs. 
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EI»0  OF  THE  SUBROCrrunE  CQRFAi- 


SW8W0UT  1  Nt  C  0  H  f  A  C  (  OUOU  ,  T  A*l,  8E  T  A  ,  M  ACM ,  Rtf  A  ,  R  AO  *  C  «0  >  I  HR  ) 

This  VERSION  OIFFCNS  IN  COHHON  AND  CALL  5>T  ATtMt^I  s 
VALUES  OF  Cl  A  NF  entered  AS  DATA 
real  HACH 

data  ci,c2,c5/  *i.iiJ,-e.eo,-o,SM/ 

FnB(UJ  ■  SQKT  (  (  U«U)  ♦  l,) 

f»*0(U)  a  l.-SQRT((M./UI*«2)/(l.»((l. /(()••?)» 

FnE(U]  *  (  C-6.*  (()••«;  .»  -  1 5  ,  »U  ®U )  -  2.)  /  {  (FNH(U>  )*•*>) 

FHF(u3»((i2o,«(u*®Fn»n2o,*(u***»))*(is<»,«(u»»un®(ioe.'*ii»u)».,a.) 

NTEhm  a  ) 

PI  *  J.l«l«* 

IF  (maCM.CT.O.BQS)  HACH  *  0.B95 
BtTA  ■  SORT  (  1.  •  <maCh«HACh)) 

SIGPB  a  CHO  •  (2,/P I ) 

AN  ■  l./(Pl*SIGPH) 

CAU  ■  BETA  •  AO 

IF  (HBA.L  t  ,0.WN)  GO  fP  uo 

Cl  8  2.0 

G2  ■  0.oN»n7 

Gi  s  o.aO 

Gw  TO  Ji 

ho  Continue 


NtAH  COHPLFHtNl A»Y  4  T  xG 

u  ■  t .  /  (¥®CAB  •  ?, , 

UP  a  l./U 

GV  a  a2.®fnu(u) 

IF  (NTfK4.NE.Jl  GO  TP  10 
G2  a  (2.  *F  NE  Cul  ♦  «.  )  't\. 

GJ  a  («ti.»(2,*f  NF(,j)/(tNH(U)»*o)  )  )/120. 
10  CONTINUE 

►AH  COMPLEMENT AHY  « I NG 


51  CONTINUE 
IF-  «  1 

IF  4  ■  0 

IF  (IF-.Nt.n  GO  TP  50 
U  ■  1,  /  <«hA*CaP»2.) 

GI  a  ♦2.*FN0(w)  ♦  Wl 

if  (ntenn.ne.jt  go  rn  jo 
G2  »  (2.»FNE(0>««.)/N.  ♦  G? 

GJ  a  (43  .- (2,  *F  nF  (li  j /(Fn«(u  )*•<»))) /I20  .  ♦  GJ 
11  CONTINUE 
JO  CUnTInuE 

OUOU  a  Cl  •  TAl  *  Gl 

IF { NT t HM . 1 0 . J)  OUOU  a  OUOu*  (C?»TAn«G2>  ♦  (C5  *  Tau  •  GJ) 
Ownu  •  OUOU  /  (  U.ovitzi' i) 

IF  (IPP.GT.O)  «KITF  (J,JOO)  Y,C  AH.r.J  ,  w2  #  G  J  ,  UtIUU,  UP  ,  0 
3oO  FORMAT  ( |M.£Fe. 5,6 ( IX, E I  1 .«) ) 


LX8TOK  OF  TSK  BUmWTTUU  TR 


to* 

4 

_# 

to 

J»4 

•to 

A 

Z) 

4* 

X 

to 

to* 

• 

ft** 

«3 

— * 

4( 

c 

4 

ft 

4 

X 

#4 

41 

(V 

tot  O 

4* 

-4 

•*4 

X 

to 

44 

c 

u 

P4 

% 

<k 

ft 

4 

w 

o 

ft 

ft 

41 

X 

to 

to* 

C» 

t 

u 

to 

ft 

to 

4 

to 

4 

U 

3 

% 

to 

a 

V 

•4 

to 

V 

7 

to* 

4 

X 

a 

V4 

4 

to* 

to 

41 

o 

to 

X 

•4 

o 

X 

u 

c 

4 

w 

C 

»■  ’ 

#-> 

4 

X 

X 

X 

% 

o 

to 

to 

4* 

a 

X 

to 

o 

W 

-* 

u 

X 

ft 

to! 

4 

4 

to 

C 

•> 

^4 

4J  ¥) 

X 

4 

>- 

•to 

4 

% 

to 

* 

X 

4 

to. 

to 

X 

c- 

to 

w 

«^» 

C- 

X 

a. 

M 

C 

to 

rvj 

»-« 

a 

-j 

to 

ft 

to 

•4 

4 

_i 

to 

to 

ft) 

4-* 

o 

to 

•0 

-X 

■4 

to* 

y 

• 

z> 

% 

_i 

4 

to 

» 

to 

O 

to 

4 

a 

to 

w 

ft 

i.‘ 

4 

4 

»- 

-I 

to 

4t 

to%  to 

k 

a 

U 

J- 

•*»  ft 

ft 

rv 

4- 

4» 

to 

ft* 

w  a 

•-• 

ft» 

to 

to 

-0 

a  ►- 

W* 

-J 

a 

4 

X 

to 

to  W 

X 

X 

a 

f- 

tr 

to 

i  a 

ftt 

ft 

— * 

4 

to 

^4 

to  UJ 

4 

IA 

u 

to 

u* 

to 

ft 

O 

»- 

«; 

» 

c 

4 

V.  ^ 

X 

to 

to* 

Kl 

to. 

*•  to 

O  *J 

u 

►- 

to 

to 

N.  _J 

4 

4 

*— . 

r> 

to 

to 

to 

m 

X 

ft 

00 

c 

A 

c  to 

C  > 

O 

A 

a 

4- 

to-  4 

2  *- 

ft* 

fti 

X 

to 

to 

ft  Z 

•  — 

41 

to 

X 

r> « 

~ 

4  to 

r  c* 

ft- 

4 

M 

tO 

c 

a  <n 

u  ft* 

ft. 

o 

•■4 

a 

ft  ftJ 

It 

—  i 
t  x 

C  N  ^ 

^  ^  C 
x  o 

C“  o  # 

K  4  /N 
2  to» 

:  i  c 

IT  *  t- 

<\l  — 


X  C*  •  •  A 

*1  (b  »>  ^ 

to  •  U  w  c  ft« 

Al  C  U  ^ 

u  A  u  ii  -r 

M  C  M 

•  w  i_3  tL  *  to 

4,  ? 

c  to.  u  -  «  '  n  o 

UhOUOU 


rt  t  u. 

r>  * 


u  u  u  u  u  u 


If  (J^T‘iT.1'0,0)  GO  to  *5 

«»TTf  (*,?<>)  »8*#'**rw#*'FrF#CriTfc''J,C0TBM,cnE^rCD 
29  F OS *•*  T  fin  ,$»,7(M9.a,?X}) 

25  CONTINUE 


LISTIHG  OF  THE  SUBROUTIHE  TIPRU 


8U1K0UTINE  TIPRLP  (XMyA,0HCM,CCSAL,VFP3,Cl.XHU,NPSI,AC,0CuTH 
A  ,RAD,CMD,TaU,3IGPN,AN.OCUTM,NSTA.NR8T,  IPR) 

C  IP  XPM  18  POSITIVE,  then  DETAIL  PRINT 
COMMON  R(30) 

REAL  MLT#ml,mino 

PI  ••  3.1414 

DtLPSI  •  2.0P1/NPST 

SUN  a  0,0 

PSI  •  0.0 

MOT  A  •  0 

MM8T  •  0 

M 100  •  (OMGN  ♦  (VPPS*COO*L) )/Cl 
IP  (MIOO.CT.O.OOS)  MRITt  (3,504)  N|«0 

304  FORMAT  (1H  .'ADVANCING  TIP  MACH  NUMBER  EXCEEDS  UN£,h( I ,00)»*F0.3) 

200  CONTINUE 

SINPSI  a  SIN (PSI ) 

MLT  a  (0M(,H*VFP8*C0SAL*5lNP5n/Cl 
IP  (MLT.LT.XMOA)  60  TO  202 
NHS  >  0 

DO  201  !•} , 30. 1 
IP  (R(I).EO.O.O)  GO  TO  202 
C 

HAV  ■  (  p  (  I  )  ♦  N(I-m/2. 

c 

ML  a  MLT  •  n»V 
IKmL.LT  ,*NOA)  GO  TO  201 
NHS  «  NR8  ♦  1 
IP  (I.PO.l)  NMITE  (3,505) 

305  FORMAT  (IN  ,»  SPAR  IS  ABOVE  ORAG  DIVERGENCE  MACH  Ny') 

IP  (ML.GT.0.M5)  ML-O.OR5 

UL  a  XMU  *  SINPSI  ♦  RAV 
GAM  »  1,  ♦  (0.2«ml*^l) 

CALL  CORFAl  (OUOU,T*U,«ETA,ML,RAV  ,NA0,CHD, IPM) 

VALUr  ■  UL*UL*GAM.OuOU«R(I)*ML*n  ,/(8tTA*«)  ) 

NS T A  a  NS T A  ♦  1 
DELS  a  MCI)  -  R(I-t ) 

Sum  a  3uh  «  (VALUE  aDElPSI *OtL») 

IP  (IPR.GT.O)  NHITF  (5,300)  PS!,R(I),ML,RAV  ,VALUt,SUM 
3CS  ( C"rii I  (r  15. a, r  i2,»,Fs,3«3(4X,E15.7)) 

IP  (R ( I ) , tO, 1,0)  GO  TU  202 

201  CONTINUE 

202  CONTINUE 

IP  (XHK.GT. 0.001)  GO  TO  214 

PS I N  a  NPSI 
SUM  ■  Sum  *  PSlN 
PSI  ■  A.2B1 
214  CONTINUE 

PSI  ■  PSI  ♦  DtLPSI 
IP  (NHS.GT.NHST)  NRST  a  NHS 
IP  (PSI. IT, s. 28)  Gn  TO  200 

OCOTR  •  SIGPB«bN  •  AC*  SUM  a  (l,/(4.*PI)) 

OCOTRS  a  OCOTH  /  <SIGPB»0N) 

*t  Turn 

DTD 


1*9 


Ml  Y2 m«>-2  -  3l*J/kN.7}  9112  PAGE.  1  _ 

C  YIP  RELIEF  CALCULATION  DRIVER 

C  _  _  _ _ 

'6ImXmSion  RR(3,30) 

COMMON  R( 3? )  _ -  - - 

REAL  *190,  LXmu 

C  IF  IPR  is  POSITIVE,  THE*  DETAIL  PRINT  _ 

IBP  FORMAT  { 4F ) 

.  .  . HI  fQHMAt.iJF) _ _ 

1B2  FORMAT  (31  ,3F, I  ) 

103  FORMAT  (4f , 2  I )  -  - 

152  FORMAT  (1M-) 

16B  FORMAT  (1*  ,1X,f?,?»3xiF8,ii?x1f10.^j2(3XiE1^.5)i1XiI>i4x,I51  __  . 

161  FORMAT  <1M-,1X, 3uxmu, 7X, 4HVFPS.8X , 7MM( 1 ,92 ) , 6X. ShOCQTR , 12X , 

A  IIHOLlCO/SIOmA j  ?* j  4wNSTAi3Xj5RNRSTA}  .  .  .  _ _ 

180  FORMAT  (INI  ,1X, 'lCaSE  IPR  vPSI'J 

181  FQRMAT  (INP.IXj 'RADIUS  Blades  OmECaR  S!G/6^ADE  'rlGMA  _  okorh/e 

2  I  A  '  > 

182  format  ( 1*0 1 2* » 'Taj  mdRAGDIV  AIRFOIL  SLOPE  SOuND^PEEO  CCSal'I 

190  FORMAT  (IR  ,1X,  13. 3X,  I?,2X,i2! 

191  FORMAT  (IN  ,lX#F6,l,lX,r4.-.*l2X,F6ll<2c2XlF9,6;i2X.,r_7J4j _ 

192  FORMAT  (IN  (F5,J,2X*F5.2»6x*r8.4,«x,F7,l,6XiF9,6} 

c 

data  (RR(1, l). 1*1. 30  l/P.?.-'. 4, 7,6,. ■.8,0,9,:’, 95. 0,97, 0,96,0, 99, 

A  0.995,1,0, 19.0,0/ 

OAT*.  (RR(2,  I  ),»  *1,30) /0, 60,0. 75,0, 90, 0,M«,  67,0,09,0  ,  VI,  KJ  ,V.U  , 

A  0.940, 0,950, 0,9550, 0.9600, 0,965  ’, 0.9700,0.9725,  __ _ 

B  P.9750, 0,9800,0. 98. '5, P.9P5T,0, 98 75, 0.9890, 0,9905 ,0,9920, 

C  ® i 9938 ,  0,9950,  0,9960,0,9970,  0,995,0,959/  ,1,0/  ..  _ 

DATA  <RR(3, I), 1*1,301/0,7,0, 75, 0,8 •’,0,85, 0,875. 0,9 00,0,925, 

A  8, 9J0,  0,960,  0.970,0,975,  "‘,98?,  0,965,  ’,990.  0  ,  V  ?5 , 1 , 0  dl4  •  .  0  /  . 

PI  •  3,1416 

1  CONTINUE  _  _ 

READ  (2,100)  RAO,  9N,  r.MGR,  CnO 

READ  (2,101)  TAU.XmUa,  Cl,  AT,  C05*L  . - 

READ  (7,102)  l  C  *5E  •  IPR,  NRSI,  L**),  Ox-M,  NhSET 

WRITE  (3,180)  -  . 

WRITE  (3,190)  lCaSE ,  IPR,  ’PR  l 

sif.eg  ■  r nq  •  _  _ _ _ 

sic  •  S 1 Gp@  •  bn’ 

x^u  ■  ux*u 

■  1,00  ♦  (i,99»oi'GR  •  nx“u/CD 

WRITE  (3.181)  • 

WRITE  (3,191)  RAO.RN.O-GR.SIGPP.SIG.CnO 

WRITE  (3,182)  - — - 

WRITE  (3,192)  TAU. xmCa,AC,C1,C0SAl 
WRITE  <  3 . 1^>2  > 

IF  (IPR,lT,0>  white  < 3,161) 

00  202  1*1,30 
202  R ( I )  *  RR(NRSET, 1 ) 

C  . . .  _ 

?0!»  CONTINUE 

VFP$  ■  XMU  •  Q"GR/COSal  .  _ 

Mive  •  (0"GR  ♦  ( vFPS#CrSAL ) ) /r i 


50 


C*u  TJPHI.F  (XH04,0HCR,C0S*l.#  VfPS.Cl.X^O.N-'Sl  ,  AC»QCOTR,WaO,C*D. 
A  UUiSUP^e^DCflTRSiNSTAi-i^T.lPRj  _  .  _ 

IP  (ir«*CT.«)  WRITE  (3,152)  _  _ 

jr  (IPR.CT,#)  wRITL  (3,161) 

WR I  Tt  (3.,16*)  *MU.VFPS,H19?#0CQT«  ,DC2TRS,;5TA,NNST _ _ _ 

X*U  •  ♦  OXNU  .  _ _ 

IP  (XMf.CT.UXHlj)  CO  TO  2?1 

IP  (P19C.6T.UMACH)  CO  TO  2/1  .  _ 

CO  TO  2*P 

2*1  CONTINUE  _  _ 

WRITE  (3,152) 

00  2*9  11*1,16,19  .  .  _ 

III  •  M  •  1* 

2*3  WRITE  (3,17*)  ( R ( I ) ,  I  •  1 1 , 1 1 1  )  _ 

'17*  FORMAT  (1*  ,  15F  7 , 4  ) 

IP  (LCASE.GT.0)  CO  TO  1  ..  .. _ 

STOP 

two 
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Change  the  equation  in  the  computer  subroutine  TR  on  page  48  from 

DELCDC  -  -  ((2.  -  MACH)/31)  *  DELMM 

to 

DELCDC  -  -  ((2.  -  MACH  *  *  2)/Bl)  *  DELMM 
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